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_ CHAPTER 1 
INTRODUCTION 


1-1. General. 
a. This series of manuals, entitled Designing 


Facilities to Resist Nuclear Weapon Effects, is 


organized as follows: 

TM 5-858-1 Facilities System Engineering 

TM 5-858-2 Weapon Effects 

T 5-858-3 Structures 

TM 5-858-4 Shock Isolation Systems - 

TM 5-858-5 Air Entrainment, Fasteners, Pene- 

tration Protection, Hydraulic-Surge 
Protective Devices, EMP Protective 
Devices , 

TM 5-858-6 Hardness Verification 

TM 5-858-7 Facility Support Systems 

TM 5-858-8 Mlustrative Examples 
A list of references pertinent to each manual is placed 
in an appendix. Additional appendixes and 
bibliographies are used, as required, for documenta- 
tion of supporting information. Pertinent 
bibliographic material is identified in the text with 
the author’s name placed in parentheses. Such 
bibliographic material is not necessary for the use of 
this manual; the name and source of publications 
related to the subject of this manual is provided for 
information purposes. 

b. The purpose of this series of manuals is to pro- 
vide guidance to engineers engaged in designing 
facilities that are required to resist nuclear weapon 
effects. It has been written for systems, structural, 
mechanical, electrical, and test engineers possessing 
state-of-the-art expertise in their respective 
disciplines, but having little knowledge of nuclear 
weapon effects on facilities. While it is applicable as 
general design guidelines to all Corps of Engineers 
specialists who participate in designing permanent 
military facilities, it has been written and organized 
on the assumption a systems-engineering group will 
coordinate design of the facilities. 

ce. Technical Manual 5-858 addresses only the 
designing of hardened facilities; other techniques to 
achieve survival capacity against nuclear weapon at- 
tacks are deception, duplication, dispersion, 
nomadization, reconstitution, and active defense. A 
facility is said to be hardened if it has been designed 
to directly resist and mitigate the weapon effects. 
Most of the hardening requirements are allocated to 
the subsidiary facilities, which house, support, and 
protect the prime mission materiel/personnel 
(PMMP). This manual is applicable to permanent 


facilities, such as those associated with weapon | 
systems, materiel stockpiles, command centers, 
manufacturing centers, and communications centers. 
d. The nuclear weapon threats considered are listed 
below. Biological, chemical, and conventional weapon 
attacks are not considered. : 
— Weapons aimed at the lone itself or at nearby: 
targets 
— A range from many, relatively small-yield 
weapons to a single super-yield weapon 
- Weapon yields from tens of kilotons to hundreds ' 
of megatons 
- Weapon delivery by aerial bombing, air-to- 
surface missile, surface-to-surface missile, or 
satellite-launched vehicle 
- Detonation (burst) of a weapon in the air, at the 
ground surface, or beneath the ground surface 
Direct-overhead bursts for a deep-buried facility 
- Near-miss bursts for a near-surface facility, pro- 
ducing peak over-pressures from tens to 
thousands of psi at the facility 
e. The designing of facilities resistant to nuclear 
weapon effects is an evolving specialty that uses a 
relatively narrow data base that incorporates both 
random and systematic uncertainties. The range of 
these uncertainties may vary from significant (order 
of 1 to 2 magnitudes) to normal (10% to 100% varia- 
tion from average values). The applicable uncertainty 
value depends on the specific weapon effect or harden- 
ing objective under consideration. Loading uncer- 
tainty is generally more significant than resistance 
uncertainty. Awareness of the appropriate uncer- 
tainty factor is essential not only for system engineer- 
ing trade-offs but in the utilization of available 
analysis or test procedures. Studies and experiments 
are being conducted to improve methodology, to better 
define random uncertainties, and to reduce systematic 
uncertainties. This manual will be revised as signifi- 
cant improvements | occur in either methodology or 
data base. 


1-2. TM 5-858-2: Weapon effects. 


This manual presents the 13 weapon effects that act 
as the facility “loads.” Each effect is described by its 
expected value, its random coefficient of variation, 
and its systematic coefficient of variation. Uncertain- 
ties in the weapon, the weapon delivery, the site, and 
other factors are considered. The effects data are 
presented both graphically and analytically. 
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CHAPTER 2 were 
NUCLEAR-WEAPON EFFECTS 


2-1. High explosive ‘and nuclear-weapon 
effects. ; 


a A nuclear explosion produces 13 potentially 
destructive effects, or environments, that must be ac- 
counted for during facility design: 

- Initial nuclear radiation 

- Electromagnetic pulse (EMP) 

- Fireball 

- Thermal radiation 

- Airblast 

- Airblast-induced ground shock 

— Crater-induced ground shock 

- Crater 

- Ejecta/debris impact 

- Ejecta/debris depth 

— Dust 

-— Firestorm 

- Residual nuclear radiation 
With the exception of the first and last, the nuclear 
radiations, these effects are also produced by a high- 
explosive (HE) explosion. Even an HE explosion can 
produce the electromagnetic pulse and the mushroom- 
shaped cloud that are commonly identified with 
nuclear explosion. The gross non-nuclear effects pro- 
duced by the two types of explosions are so similar 
that the explosive energy of a nuclear explosion is also 
expressed in terms of an equivalent TNT-charge 
weight. 

b. This equivalent TNT-charge weight is referred 
to as the explosive energy release, or yield, W. A facil- 
ity designed to withstand the effects produced by an 
HE explosion would be quite similar to a facility 
designed to withstand the effects produced by a 
nuclear explosion of equivalent yield. This statement 
is true for peak overpressure larger than a few tens 
of psi, which are the hardness levels addressed in this 
volume. 

c. Since the protective measures required for the 
non-nuclear effects also provide, in general, sufficient 
protection against nuclear radiation, the protective 
measures required for the two facilities would be of 
the same kind—but would differ drastically in degree 
when the effects are produced by weapons of war. HE 
weapons (e.g., conventional bombs) have yields of only 
several tons or less, whereas state-of-the-art nuclear 
weapons range in yield from a few kilotons to a few 
tens of megatons. The potential destructiveness of a 
weapon at a given distance varies between W!/3 and 
W, depending on the particular weapon effect under 
consideration. Thus, a small state-of-the-art nuclear 
weapon is 10 to 1000 times more destructive than the 
largest HE weapon, and a large nuclear weapon, 1000 


to 1,000,000 times more destructive. It is this much 
larger magnitude of potential destructiveness that 
creates disparity between hardening against nuclear 
weapons and hardening against conventional HE 
weapons. 


2-2. Nuclear weapons. 


Nuclear weapons include atomic and thermo-nuclear 
weapons. The explosive energy of an atomic weapon 
results from the fission of uranium and plutonium; 
the yields of atomic weapons are tens of kilotons or 
less. A thermo-nuclear weapon employs both fission 
and fusion reactions to obtain a very large yield. An 
ordinary nuclear weapon of either basic type can be 
modified to produce an “‘enhanced’’ weapon, e.g., the 
output of the weapon can be channeled to produce 
relatively large amounts of nuclear radiation, or if the 
weapon is so designed, it may enter the earth prior 
to detonation resulting in a buried or partially buried 
detonation. 


2-3. Burst elevations. 


a. The burst elevation of a nuclear weapon can be 
classifed as follows: 
100,000 ft c High-altitude 


burst 
3500 ft/Mt®4 =< Air burst < 100,000 ft 
0 Low-altitude < 3500 ft/Mt0-4 
burst 
Contact burst = 0 
1000 ft/Mt0.25 < ze aaa Z 0 
urst 
Deep-buried — < _1 000 stamt0-25 
burst 


Note the following characteristics. Air burst: Fireball 


does not touch the ground. Contact burst: Center of 


gravity of weapon is slightly above the surface; some 
penetration required with true surface burst in which 
c.g. is at ground surface. Deep-buried burst: Com- 
pletely contained and does not vent to the 
atmosphere. , 

b. State-of-the-art weapons cannot penetrate to suf- 
ficient depths to result in a deep-buried burst, 1.e., a 
contained burst. High-altitude and air bursts general- 
ly are not effective burst elevations for the higher 
hardness levels considered in this volume. However, 
there may be threat scenarios that would require the 
consideration of these burst elevations; their effects 
are selectively included in this volume. Low-altitude 
bursts, contact bursts, and shallow-buried bursts pose 
the greatest threat. They are collectively referred to 
in this volume as near-surface bursts. 
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2-4. Weapon effect, a random variable. 


a. Parameters. The specification of a weapon effect 
for use during the designing of a facility requires that 
the following basic parameters be quantified: 

— Weapon design 

— Weapon yield 

- Burst elevation 

- Slant range | 

- Weather and atmospheric conditions 

-— Geologic factors : 

Because of the Ge aie with these basic 
parameters, it is necessary to treat them, and thus 
the weapon effects, as random variables. 

b. Definition of a random variable. A random 
variable-X is a numerical variable, the specific value 
of which cannot be predicted with certainty. Every 
specific value of the random. variable X, denoted by 
x, represents a possible outcome of the random 
phenomenon. under consideration. 

¢. Density function. The probability density func- 
tion of X, denoted by fy. (x), describes the relative 


likelihood of the various outcomes of X. It is defined — 


such that . 
fiss00 dx = 1.0 
oO . : : 


where f(x) dx represents the likelihood of X being 
between x and x + dx. This relative likelihood or 
relative frequency of occurrence is known as the pro- 
bability of x < KX <x + dx,ie, 


(2-1) 


Px < X < x + dx) = fx(x) dx (2-2) 


Also, note the probability that x, < X < x» is 


Pex, <X< Xo) - {* fx(x) dx (2~3) 
x 


1 


d. Mean and variability ofarandom variable. The 
mean (the average, or expected) value of the random 
variable X is 


By =f ae fy(x) dx 


A useful measure of the randomness, or variability, 
of the random variable X is its coefficient of variation: 


(2-4) 


dy = Ox/py (2-5) 
where ox is the standard deviation of X and is equal 
to ee 


oa [J 7 oz Hy)? fy dx 12 
oo ‘ 


(2-6) 


e. Evaluation of probabilities. In general, obtain the 
probability of x, < X < x» from equation 2-3. For a 
number of special density functions, these pro- 
babilities are tabulated in reference books. (see, for 
example, Abramowitz, 1972). As an example, for a 
Gaussian or normal density function 


1 
of or 


Letting z denote the scribe: af standard deyienone 
above mean, one obtains 


fy(x) = exp | -0.5(x - Hs. 08 (2-7) 


zy + 20x 
P(z) = P(X uy + Zoy) = a exp 
: a Oy Yan . : 
é 2 
a u 
0.5 ( Hgts exp| | du 
[- Ree [fe ab 


(2-8) 


where transformation u = (x — py) /ox is used. Values 
of this integral for a number of values of z are given 
below. 


z Piz) 

0 0.500 
0.1 0.540 
0.25 0.599 
0.50 0.691 
0.75 0.773 
1.00 0.841 
2.50 © 0.994 
5.00 0.999 


The above table may be used to evaluate other pro- 
babilities for the Gausian distribution; for example: 
P(-z) = 1 - Piz) 
Ply + 2)0x <X< py + Z Oy) = Plzo) - Plz,) 

f Accounting for uncertainties: A single variable. 


The uncertainty associated with a random variable 
is due to two factors: (1) the natural or inherent ran- 


domness of the random variable itself, and (2) errors — 


which result from the inability to accurately estimate 
the mean, the standard deviation, or the density func- 
tion of the random variable. The latter may result, 
for example, from lack of sufficient-data. To account 
for both uncertainties the random variable X is writ- 
ten as 


X =N,X (2-9) 
where X is the random variable as represented by the 
available data and has a mean x and a coefficient of 
variation dy which are estimated from the data. Ny 
is a random correction factor which is included to 
model the estimation errors. Its mean value Ny 
represents the bias in data, whereas its coefficient of 
variation Ay is a measure of the uncertainty 
resulting from estimation errors. As a first approx- 
imation Ay is considered as the error in estimation 
of the mean value x. The mean value of the random 
variable X is evaluated from 


vy = Nyx (2-10) 


However, if an unbiased set of data is used to evaluate 
x, then Ny = 1.0 and 
i.e., for an unbiased set of data, x is the best estimate 


of the actual mean. The total coefficient of variation 
of X, designated by y, is obtained from 


ae 
Qy = V Oo. AS 
Here, dy models the natural randomness of X, 
whereas Ay represents the uncertainty arising from 
errors in estimation. As an example, when a set of 
n unbiased samples are available, it is usually as- 
sumed that 


(2-12) 


g. Accounting for uncertainties; Multiple variables. 
If X itself is a function of several random variables, 
De TY Voy Ngee eige Vad ve WD (2-14) 
then using first-order linear approximations (Ang- 
Cornell, 1974) 


K = f(¥), Yo, ¥s, Yas , Yn) (2-15) 
and 
N y? Z 
‘ of 
Q2 = Q2 + = Q2 
x f 2 x2\ 0Y, 4 Yy 
| af af ase 
7 “>2 | OY. OY Onm**Y,,*Y 
n#¥m mn ° ‘ - 
(2~16) 
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where s 
Q- = Total coefficient of variation 
associated with the functional form of 
xX ‘5 
= 2 A2 
ds As 
Qy = Total coefficient of variation 
associated with Y,, 
2 2 
Oy. + a 
Cnm = Correlation coefficient of Y, and Y,, 


~-1< Cam < <1 
The subscript “o” in equation 2-16 denotes that 
6f/AY,, is to be evaluated at the mean values of the 
variables. Implicit in equations 2-15 and 2-16 are the 
assumptions that 2% << 1 and the nonlinearity in 
f near the mean valubs i is not large. In this manual, 
both symbols 2 or 22 are called “uncertainty” to - 
avoid the repetitive use of the term “uncertainty 
squared.” 

h. Distributions, If a particular distribution is de- 
fined, the uncertainty can be calculated in terms of 
the upper and lower bounds assigned to the data. For 
example, suppose that data are presented within 
bounds that are said to include 95% of the data. If the 
data are assumed to have a Gaussian distribution, the 
bounds correspond to the + two-standard deviation 
limits (see e¢ above). Procedures for a number of 
specific distributions are described in the following 
paragraphs. 

(1) Normal: The coefficient of variation of a nor- 
mal distribution is 


qo 1l;R-1 
xX n\|R+1 
where R = L,/L,, and L, and Ly are the lower and 
upper bounds corresponding to n standard deviations 

on either side of the mean. 


(2) Log-normal: The coefficient of variation of a 
log-normal distribution is 


7 (InR? ue 
é>[ oe [SR] ] 


where R = k?, k = X/L, or k = Lo/x, and Ly and Ly 
are the lower and upper bounds of the distribution. 
The quantity n is that number of standard deviations 
of the bounds from that log mean. 

(3) Uniform: The coefficient of variation of a 
uniform distribution is 


(2-17) 


(2-18) 


(2-19) 
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where L, and Ls are the lower and upper bounds of 


the data. 
(4) Log-Uniform: The coefficient of variation of a 
log-uniform distribution is 


& = (g>2)\ io 


where L, and L, are the lower and upper bounds of 


the data. 
i. Uncertainty in nuclear weapon effects. Uncertain. 


V2 


(2-20) 


ties in nuclear weapon effects arise from yield, 
delivery accuracy, weapon capability (air bursts, 
penetrating bursts), as well as from imperfect 
knowledge of specific effects for known yield, range, 
etc. In the following chapters an attempt is made to 
identify the level of uncertainty of our knowledge of 
weapon effects. The uncertainties relative to weapon 
yield, range, and burst conditions are defined, but 
since they are considered a part of threat definition 
they are not quantified; however, provision is made 
for incorporating uncertainties in these parameters 
from threat definition applicable to specific problems. 
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CHAPTER 3 . 
INITIAL NUCLEAR RADIATION 


3-1. Gross phenomenology. 


a. A nuclear explosion is accompanied by the emis- 
sion of nuclear radiations consisting of gamma rays, 
neutrons, beta particles, and alpha particles. The 
range of beta and alpha particles is so short that they 
can be neglected in the estimate of the initial radia- 
tion. The initial nuclear radiation is considered to be 
composed of neutron and gamma ray radiation from 
both initial effects (t < 1 minute) and secondary in- 
teractions. Because induced radioactivity (a threat to 
personnel for many hours) begins during t < 1 min, 
it is included in this chapter. 

b. The basic variables affecting the destructiveness 
of initial nuclear radiation are: 

- Weapon design 

— Weapon yield 

- Burst elevation 

~ Slant range 

— Atmospheric conditions 
The destructiveness of the initial nuclear radiation 
will be enhanced by any of the following: Larger 
weapon yield, higher burst, smaller slant range. 

c. The primary direct effect of nuclear radiation is 
an antipersonnel effect. Direct nuclear radiation ef- 
fects on equipment and materials are less significant, 
but certain detector materials and electronic com- 
ponents can be damaged. The designer must either 
design for these effects (by shielding, for example) or 
he must preclude them (by burial, for example). The 
two significant aspects of initial nuclear radiation, 
neutron radiation and gamma ray radiation, are con- 
sidered separately. 


3-2. Neutron radiation. 


a, Initial. The weapon type will have an influence on 
the energy spectrum of the source radiation. The loca- 
tion of the weapon source (and the target) will also 
affect the spectrum. (Specific situations relating these 
parameters are too numerous and specialized for 
presentation here; however, such information can be 
found in the report prepared by SAI, 1971). Only in- 
itial neutron radiation effects (occurring within 1 
minute after detonation) are considered in this 
chapter; residual neutron effects (after 1 minute) are 
discussed in chapter 15. 

b. Equipment: The vulnerabilities of electronic 
equipment are expressed in terms of fluence spectra 
(i.e., incident neutrons/in.? or neutrons/cm? at various 
energy levels expressed in MeV). Typically, fluence 
spectra are presented in charts of the kind shown in 
figure 3-1, which are based on the weapon (source) 


spectrum and the attenuation characterization of the 
standard atmosphere. 

c. Personnel. Personnel injury resulting from radia- 
tion exposure is complex and numerous (Glasstone, 
1977). If the type of weapon treat is known and if the 
slant range is less than 5000 yd (4572 m), it is possi- 
ble to quickly estimate the neutron does to person- 
nel regardless of the spectral content of the pulse. For 
low-altitude atmosphere bursts, where personnel are 
located on or near the earth’s surface, the neutron 
dose absorbed by human tissue is (DNA, 1972). 


WFYH, 
Dy = ——— rad (J/kg) . (3-1) 
R2 
where 
W = Weapon yield, kt 
R= Slant range, yd (m) 
Fy = Normalized neutron dose, yd? 
rad/kt[m? (J/kg)/kt] 
H, = Burst height correction factor 


(1) The mean value of equation 3-1 is calculated 
by referring to figures 3-2 through 3-5 (which show 
the relationship between Fy and R), and to figure 3-6 
(which shows the relationship between H, and the 
height of burst) for the standard division (DNA, 1972) 
of weapon types, I through VIII. 

(2) The uncertainty associated with equation 3-1 


2 = 2 2 2 2 2 3-2) 
Q QF + Qi + ME + QE + 4 Bp (3-2) 


where 8, is the uncertainty of the functional form of 
equation 3-1, Qy is the uncertainty of the yield, etc. 
Consider that 2; = 0.2, Qp, = 0.30; using equation 
2-19 and the assumption that figures 3-2 through: 
3-5 represent mean values of data uniformly 
distributed with upper and lower bounds within 
+25% of the mean (from DNA, 1972) and Qy, = 0.3. 
Thus, equation 3-2 reduces to 


ae 2 2 3-3 
ce 0.22 + 92, + 4 92 (3-3) 


3-3. Gamma radiation. 


a. Prompt and delayed. Gamma radiation can be 
divided into two components: Prompt (within 1 or 2 
shakes, a shake being equal to 10-° sec), and delayed 
(up to 1 min). Electronic equipment is sensitive to the 
energy spectrum of the radiation—only prompt dosage 
is calculated; personnel are influenced primarily by 
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Figure 3-1. Typical Neutron Fluence Inetdent on a Receiver 
Located on or near the Surface of the Earth (DNA, 
1972) 
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Figure 3-2, Neutron Dose Fy, as a Function of Slant Range from a 1-kt Surface Burst, Weapon Types I through IV. Short Ranges (DNA, 
1972) 
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Figure 3-3. Neutron Dose Fy as a Function of Slant Range from a 1-kt Surface Burst, Weapon Types | through IV, Long Ranges (DNA, 
1972) 
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Figure 3-4. Neutron Dose Fy as a Function of Slant Range from a 1-kt Surface Burst, Weapon Types V through VHI, Short Ranges (DNA, 
1972) 
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Figure 3-5. Neutron Dose Fy as a Function of Slant Range from a 1-kt Surface Burst, Weapon Types V through VIII, Long Ranges (DNA, 
1972) 
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Figure 3-6. Burst Height Adjustment Factors H, for Neutrons and Secondary Gamma Rays (DNA. 1972) 


the total dose absorbed—both prompt and delayed ef- 
fects are important because the effect is cumulative. 

b. Prompt gamma ray dose on equipment. Prompt 
gamma radiation that effects electronic equipment 
can be expressed by an equation of the form 


K,W 
7 es exp (KyeR) rad (J/kg) (3-4) 
where 
. K,, Ky = Constants 
= Yield in kt 
R = Slant range in cm 


e = Air density in gm/em? 
Other sources of gamma rays include inelastic scat- 
tering of neutrons and the decay of fission products. 
These sources depend primarily on weapon design 
(SAI, 1971). Equation 3-4 will generally apply to 
materials with atomic number 20 or less. 


(1) The magnitude of the constants in equation 
3-4 should be determined from the upper bound 
estimates of the prompt gamma ray flux equation 
(o,,,). Moreover, the source energy density MeV/cm? 
should be converted to absorbed erergy expressed in 
rays (by the conversion factor presented in the 
reference), and the result multiplied by 2 shakes to 
obtain the prompt gamma ray dose. 

(2) The uncertainty associated with equation 3-4 
is expressed by 


2 =~ 924 92 2 2 p2R2 Q2 2 02 
Se Et 2 + Sw t Kee R 2° + (RRg@ - 2) QF 
(3-5) 


where ©, is the uncertainty in the form of equation 
3-4, Qy is the uncertainty in the yield, etc. Assume 
that Q, = 0.2, Q, = 0.05, and that Q,. = 0.3 (using 
equation 2-19 and the assumption (DNA, 1972) that 
the bounds are within + 33% of the mean). 
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c. Secondary gamma ray dose on personnel. In ad- where 
dition to the neutron dose described in paragraph 3-2 
c, secondary gamma rays and fission-product gamma 


rays will be significant to personnel. If the type of W = Weapon yield, kt 

weapon is known and if the slant range is less than Be Normalized secondary gamma ray 
5000 yd (4572 m), it is possible to estimate the secon- dose, yd? rad/kt (m2(J/kg)/kt) 

dary gamma ray dose. For a low-altitude burst in a H, = Detonation height correction factor 
standard atmosphere where the personnel are located R = Slant range, yd (m) 


on or near the earth’s surface, the secondary gamma 


Pay Goee (1) The mean values of equation 3-6 are 


WF_H | calculated by referring. to figures 3-6 through 3-10 
Ds ee) oad (J/kg) (3-6) for the standard division of weapon types (DNA, 


R2 1972), I through VIII. 
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Figure 3-7. Secondary Gamma Ray Dose F,, as a Function of Slant Range from a 1-kt Surface Burst, Weapon Types | through IV, Short 
Ranges (DNA, 1972) 


3-8 


Figure 3-8. 


Fys AIR GROUND SECONDARY GAMMA RAY DOSE, ya? rad/kt 


SLANT RANGE, m 


T™ 5-§ 


2000 2500 3000 3500 4000. 4500 
Peel: Nee any) 


RANE 


9! 


r=) ° 
Baal wh 
T Trt Pe TTT 
S ° 
i) a : 
AIR GROUND SECONDARY GAMMA RAY DOSE. m@(J/ka)/kt 


EEE aged Coren 


107) 


2000 . 2500 3000 3500 4000 4500 5000 


R SLANT RANGE, yd 


Secondary Gamma Ray Dose se as a Function of Slant Range from a 1-ki Surface B 
Ranges (DNA, 1972) 


urst, Weapon Types I through IV, 


TM 5-858-2 


SLANT RANGE, m 
0 500 _ 1000 1500 2000 2500 


7 
£ 10? 107 
=S = x s 
@ a = 
L m 
N dé 
Bd r S 
a 6™ 
w ..8 10 E 
2 10 Eo ss 
a is 7) 
> = oO 
= L. a 
L > 
: = 
< 5 < 
& 40/ 1° = 
> pr 
a = 
3 a 
S 10 E10 Ye 
= - fas) 
2 E z 
3 f r= 
wz - a 
_ L o 
< a 
n ae * 
uP 195 10 ? 
C 
= 
tol L-10 


) 500: 1000 1500. 2000 2500 3000 
R SLANT RANGE, yd 


Figure 3-9. Secondary Gamma Ray Dose F,, as a Function of Slant Range from a 1-kt Surface Burst, Weapon Types V through VII, 
Short Ranges (DNA, 1972) 
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Figure 3-10. Secondary Gamma Ray Dose F,, as a Function of Slant Range from a 1-kt Surface Burst, Weapon Types V through VIII, 
Long Ranges (DNA, 1972} . 
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(2) The uncertainty associated with equation 3-6 
is 


+ Q2 


2 = 2 2 2 
Q2 OP Ok Sd OF ae 


+ 4Q2 3-7 
2 (3-7) 


ys 


where 2, is the uncertainty in the form of the equa- 
tion, Qy, is the uncertainty of the yield, etc. Consider 
that Q- = 0.2, Qp , 7 08 and using equation 2-19 
with the assumption that figures 3-7 through 3-10 
represent mean values of data uniformly distributed 
with bounds of +25% (DNA, 1972), and 2H, = 0.3. 
Thus, equation 3-7 reduces to 


25 = 0.22 + 92, + 402 (3-8) 
d. Fission-product gamma ray dose on personnel. If 
the weapon yield, height of burst, slant range, and 
air density are known, the fission-product gamma ray 
dose absorbed by human tissue can be estimated. The 
fission-product gamma ray dose D,, is 


where 
W, = Fission yield (kt) 
Ee = Normalized fission-product gamma 


ray dose, rad/kt [(J/kg/kt)] 
Hp = Slant-range-dependent adjustment 


factor 
Hwy = Yield-dependent adjustment factor 
E = Hydrodynamic enhancement factor 


(1) The mean value of equation 3-9 can be 
calculated by referring to figures 3-11 through 3-20. 
In these figures, @ = o/e9, where @ and @, (= 1.2250 
x 10-3 gm/em3) are, respectively, the mean actual 
and ambient standard air densities (15.5°C and 14.7 
psia). 

(2) The uncertainty of equation 3-9 is 
Qh = Q2+ i + i oF + Qi + Q2 (3-10) 
where &; is the uncertainty of equation 3-9, Sv, is 
the uncertainty of the yield, etc. Using equation 2-19 
and the assumption that figures 3-11 through 3-13 
represent mean values of data that are uniformly 
distributed and whose bounds are within +25% of the 
mean values (DNA, 1972), and assuming other values 
of uncertainty, 


Rp = 0.2, QP = 0.3, and Sip = QHy = Qp = 0.3; 
then equation 3-10 reduces to 


—_— 2 
OF, = 04 + Qi (3-11) 


Dyf 


3-12 


e. Total personnel dose. The total dose to personnel 
is the sum of the neutron dose, the secondary gamma 
ray dose, and the fission-preduct gamma ray dose. 

(1) The total mean dose to personnel is the sum 
of equations 3-1, 3-6, and 3-9. 

(2) The uncertainty Qe of the total dose to person- 
ne] is the uncertainty of the sum of equations 3-1, 
3-6, and 3-9. Thusly, 


H2w? 
2R4 


i 2 
27 = 0.09 + Qt 


2 2 2 
fo09 F2 + (0.09 + 402) (Fy + Fy.) | 


2 24 2492 Pr2 
_ 036 WEFEHRHE 


2 
Dy 


(3-12a) 


except for megaton weapons where R < 1000 yd 
(914 m), then 


H2w2 
Q2 = 0.17 + Of + 1 
D2R4 


2 _ 402 2 
| os F2 + (0.09 + 422) (Fy + F,,) 


202 FI2q2 K2 
0.36 WeF HR AWE 
D2 


+ {38-12b) 


f Neutron-induced gamma activity in soils, residual 
effect on personnel. Neutron-induced soil activity must 
be considered for a scaled burst height below 75 
ft/kt®35 (23 m/kt9-35) but above the ground. Soils that 
contain sodium, manganese, and aluminum will con- 
tribute most to the induced radioactivity, For times 
after burst of less than one-half hour, aluminum is 
most important; between one-half hour and five 
hours, manganese is most important; at later times, 
sodium becomes increasingly more important. 


(1) Figure 3-21 presents upper and lower bounds ~ 


on the neutron-induced gamma dose rate at a 
reference time one hour after burst. At any given 
slant range, assume that the mean is the logarithmic 
average of the upper and lower bounds. For a time 
other than one hour after burst, this mean must be 
multiplied by the factors shown in figure 3-22 for the 
four types of soil presented in table 3-1. Additionally, 
all data must be multiplied by the factors shown in 
table 3-2 to obtain a final estimate of the mean dose 
rate. 
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Figure 3-11. Fission Product Gamma Ray Dose F,;as a Function of Slant Range from a 1-kt (Fission Yield) Surface Burst, Short Ranges 
(DNA, 1972) 
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Figure 3-12. Fission Product Gamma Ray Dose F,,as a Function of Slant Range from a 1-kt (Fission Yield) Surface Burst, Intermediate 
Ranges (DNA, 1972) j 


10 


TM 5-858-2 


SLANT RANGE, m 
3000 3500 4000 4500 5000 I se ee 


INNNOSE 
IANANN \ 


ad 
= re 
a=) = 
8 0.8 5 
‘ om 
Z < 
s | re 
2 l | ~ 14 
% Pa 
4 4978 : 17? 3 
= , : 0.9 ©. 
[aa ~< 
< C e 
3 : 
a = 
a -¢ 
vas 10 7 : 10 - = 
1.0 

a 

pl 

4 

108 : 107! 
1d 


-9 Me | | ert 
10 000. 3500 “h000 ~—«4500 5000 5500 6000 
- R SLANT RANGE, yd 


Figure 3-13, Fission Product Gamma Ray Dose Fas a Funetion of Slant Range from a 1-kt (Fission Yield) Surface Burst, Long Ranges 
(DNA, 1972) 
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Figure 3-14. Range-Dependent Burst Height Adjustment Factors Hp for Fission Product Gamma Rays (DNA, 1972) 
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Figure 3-15. Yield-Dependent Burst Height Adjustment Factors Hy for Fission Product Gamma Rays (DNA, 1972) 
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Figure 3-16. Fission-Product Gamma Ray Hydrodynamic Enhancement Factors E as a Function of Slant Range for Relative Air Density 
of 1.1 (DNA, 1972) 
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Figure 3-17. Fisston-Product Gamma Ray Hydrodynamic Enhancement Factors E as a Function of Slant Range for Relative Air Density 
of 1.0 (DNA, 1972) 
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Figure 3-18, Fisston-Product Gamma Ray Hydrodynamic Enhancement Factors E as a Function of Slant Range for Relative Air Density 
of 0.9 (DNA, 1972) 
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Figure 3-19, Fission-Product Gamma Ray Hydrodynamic Enhancement Factors E as a Function of Slant Range for Relative Air Density 


_ of 0.8 (DNA, 1972) 
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Figure 3-20. Fission-Product Gamma Ray Hydredynamic Enhancement Factors E as a Function of Slant Range for Relative Air Density 
ae of 0.7 (DNA, 1972) : 
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Figure 3-21. Neutron-Induced Gamma Dose Rate D,; as a Function of Slant Range at a Reference Time of 1 Hour After Burst (DNA, 1972) 
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Figure 3-22. Decay Factors for Neutron-Induced Gamma Activity, D,; (DNA, 1972) 


3-24 


t 
_ 


DECAY FACTOR 


T 5-858-2 


Table 3-1. Soil Types by Composition (DNA, 1972) 


Percentage of Soil Type (by weight} 


Type I 
(Liberia, 
Element Africa) 


Sodium 


Manganese ~ 


Aluminum 
Tron 
Silicon 
Titanium 
Calcium 


Potassium 


' Hydrogen 
Boron 
Nitrogen 
Sulfur 
Magnesium 
Chromium 


Phosphorous 


Carbon 


Oxygen 


Multiplying 
Factor 


(2) Ina particular situation, it may be necessary 
to construct a representative decay factor curve based 
on the chemical composition of the soil as compared 
to the data presented in table 3-1. For example, 
assume that a particular site has properties in which 
the percentage of aluminum most nearly matches 
that of Soil I but that manganese and sodium are most 
accurately represented by Soils II and ITI, respec- 
tively. Construct a decay factor curve from figure 
3-22 which approximates Soil I for times up to one- 


Type II 
(Nevada 
desert) 


Type IV 
(beach sand, 
Pensacola, 
Florida) 


Type III 
(lava, clay, 
Hawaii) 


half hour, Soil II from one-half hour to five hours, and 
sodium after four hours. 

(3) Obtain the total gamma ray dose experienced 
by personnel from figures 3-23 through 3-26, which 
are applied to the one-hour dose rate presented in 
figure 3-21. The factor obtained from the appropriate 
soil type in figures 3-23 through 3-26 depends on the 
time of entry into a contaminated area and the 
amount of time spent in that area. These factors are 
multiplied by the appropriate dose rate in figure 3-21. 

(4) As an example, for the particular soil type 
discussed in (2) above, assume that the area is to be 
entered after two hours and to be occupied for seven 
hours. Enter figure 3-24 (for Soil UD) for an entry time 
of two hours and an occupancy time of three hours 
(period during which manganese contribution is most 
important). Read the multiplying factor. Then enter 
figure 3-25 (for Soil III) for an entry time of five hours 
and occupancy time of two hours (period during which 
sodium is most important) and read the multiplying 
factor. Add the factors and multiply this number by 
the mean value of data in figure 3-21. 
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Figure 3-23. Total Radiation Dose Received frad/kt) in an area contaminated by Neutron-Induced Gamma Activity, Sotl Type I(DNA, 1972) 
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Figure 3-24. Total Radiation Dose Received frad/k) in an area contaminated by Neutron-Induced Gamma Activity, Soil Type H (DNA, 1972) 
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Figure 3-25. Total Radiation Dose Received frad/kt) in an area contaminated by Neutron-Induced Gamma Activity, Soil Type II (DNA, 1972) 


3-28 


10 


; INFINITE TIME 


MULTIPLYING FACTOR, APPLY TO 1-HR DOSE RATE 
° 
fo) 


10° 


ENTRY TIME, hours after detonation 


Figure 3-26. Total Radwation Dose Received (rad‘kts in an area contaminated by Neutron-Induced Gamma Activity, Soul Type 1V t 
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(5) As an estimate of uncertainty of the neutron- 
induced gamma activity in soils Qp,, use the 
equation. - : 


Dy & © eae R Re (3-13) 


where Q, is the uncertainty of total dose and Qp is 


the uncertainty in slant range (yd). Assuming that 


the dose rate is log-uniformly distributed in figure 
3-21, calculate 2; from equation 2-19 where L, and 
Ly are the lower and upper bounds at:any slant range 
in figure 3-21. . 
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ee ELECTROMAGNETIC PULSE (EMP) a neha a” 


4-1. :Gross phenomenology. 

a. Time-varying electric and magnetic fields are pro- 
duced by nuclear detonations. Their frequency spec- 
trum varies from direct current to about 10 GHz. The 
gamma radiation produced by the nuclear detonation 
impinges on air molecules and induces a Compton 
current via ‘ionization of the air. As the Compton cur- 
rent expands away from the source, separation of the 
ions occurs and produces a radial electric field in the 
entire deposition region. Because of the atmospheric 
density gradient, the earth’s magnetic field, and the 
weapon configuration itself, a magnetic field is also 
created, resulting in the radiation of an elec- 


tromagnetic wave for a medium-altitude (where the — 


deposition region does not touch the earth) or a high- 
altitude burst. In a surface burst, the relatively large 
conductivity of the earth induces electric and 
magnetic fields (as shown in figure 4-1) whose inten- 
sities are much larger than for an air burst. 

b. The basic variables contributing to the EMP for 
exposed, flush, and shallow-buried facilities are 

— Weapon design 

— Weapon yield 

- Burst elevation 

— Atmospheric conductivity 

- Ground conductivity 


DEPOSITION 
REGION BOUNDARY 
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m 
ul 


m 
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Figure 4-1. 


c. Hardened ground-based or shallow-buried 
systems are subject to EMP caused by high-altitude 
bursts or by contact (or penetrating) bursts. Each 
system must be investigated to determine which burst 
condition predominates. Hardened ground-based or 
very shallow-buried systems will be susceptible to 
EMP caused by contact (or penetrating) bursts only 
if within the deposition region. The deposition region 
is that region in space where EMP fields are 
generated. 


4-2. Surface burst. 


a. Ina surface burst, the conducting ground provides 
an effective return path for electrons traveling out- 
ward in the air from the burst point. These current 
loops produce very large azimuthal magnetic fields 
in the deposition region, especially near the ground. 
Large electric fields are also produced. The EMP for 
a surface burst will be one to two orders of magnitude 
greater than for air bursts. 

b, The magnitude and shape-of the electric and 
magnetic fields in the deposition region cannot be 
described by simple functional relationships because 
of the complex nature of the underlying variables. For 
maximum accuracy, the currents induced in a system 
should be calculated with the recognition that the 
system, itself, influences the characteristics of the 


MAGNETIC FIELD 
RADIAL ELECTRIC FIELD 


TRANSVERSE ELECTRIC FIELD 


Surface Burst EMP 
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EMP environment. Exact solutions for individual 
cases can be solved by using the computer programs 
described in the DNA EMP Handbook (GE-TEMPO, 
1971) and in the NWE Computer Code Directory (GE- 
TEMPO, 1974). For approximations, the relative 
changes in peak values of both the radial and polar 
electric fields and the azimuthal magnetic fields 
created in the deposition region are presented in 


By PEAK MAGNETIC FIELD, Wb/m~, LOG SCALE 


figures 4-2 through 4-4 for two weapon yields and 
two ground conductivities. The time histories for 
these same fields are presented in figures 4-5 through 
4-7. Their corresponding frequency spectra are 
presented in figures 4-8 through 4-10. These data 
should be used for guidance only. Situations involv- 
ing other weapon yields or conductivities will require 
the execution of the appropriate computer program. 


aa OVERPRESSURE, psi LOG SCALE 


Figure 4-2, Peak Magnetic Field B, Vs. Ouerpressure at the Air/Ground Interface for Varying Ground Conductivities and Yields (DNA, 


1972) 
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Figure 4-3. Peak Radial Electric Field Ep Vs. Overpressure at the 
Air/Ground Interface for Varying Ground Conduc- 
tivities and Yields (DNA, 1972) 
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Figure 4-4. Peak Transverse Electric Field E, Vs. Ouerpressure 
at the Air/Ground Interface for Varying Ground Con- 
ductivities and Yields (DNA; 1972) 
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Figure 4-5. B,-Time Waveform at the Air/Ground Interface for Several Overpressure Levels (SAL 1971) 
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ER RADIAL ELECTRIC FIELD, V/m, LOG SCALE 
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Figure 4-6. E,-Time Waveform at the Air/Ground Interface for Several Ouverpressure Levels (SAI, 1971) Sy 


cy 


Eo TRANSVERSE ELECTRIC FIELD, V/m, LOG SCALE 


t TIME, s, LOG SCALE 
Figure 4-7. E,Time Waveform at the Air/Ground Interface for Several Overpressure Levels (SAI, 1971) 
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Figure 4-8 Fourier Amplitude of B, Waveform at the Air/Ground Interface for Several Overpressure Levels (DNA, 1972) 
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Figure 4-9. Fourier Amplitude of E, Waveform at the Air/Ground Interface for Several Overpressure Levels (DNA, 1972) 
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[E, (f) | AMPLITUDE, V/m/Hz, LOG SCALE 


f FREQUENCY, Hz, LOG SCALE 
Figure 4-10. Fourier Amplitude of F, Waveform at the Air/Ground Interface for Several Overpressure Levels (DNA, 1972) 


c. Facilities buried in the earth will experience a 
decrease in the amplitude of the EMP fields that oc- 
cur at the air/earth interface. Facilities elevated 
several meters above the surface will experience a 
general (but not universal) attenuation of the EMP 
fields. Figures 4-2 through 4-10 illustrate 
characteristic curves. The uncertainty is derived as 
described in d below. 

d. The uncertainties in predictions for surface-burst 
EMP depend on many variables; but for the threats 
of interest in this volume, the prime factors are over- 
pressure and weapon yield, frequency, and calcula- 
tional procedures. These uncertainties are sometimes 
referred to as tactical and calculational errors (AFWL 
EMP Handbook, 1972), which correspond to the terms 
“random” and “ignorance” uncertainties as used in 
this manual. Typical uncertainty bands for certain 
conditions are presented in figures 4-8 through 4-10. 
Typical examples are presented in figures 4-11 
through 4-13 for the uncertainty bands on the 
Fourier amplitudes of the By, Ep, and E, fields at the 
Pyo. range at a yield of W, and conductivity of 10-2 
mho/m. (Similar curves are available in the AFWL 
EMP Handbook for other overpressure regimes and 
yields.) These data may be used directly to assess the 
uncertainties of EMP environments. 


4-10 


e. The method used for determining uncertainty 
throughout this volume is embodied in equation 2-16. 
Because of the lack of useable information, it is 
customary to neglect the second term (relating cor- 
relation between uncertainties). Thus, 


where K = g(¥), Vo, . 


.. ¥,) has the same meaning as 
in equation 2-15. Wherever possible in this manual, 
the function g is defined in closed form so that the 


derivatives ag can be directly evaluated for each 


independent variable (such as P,,, W,, Hz, Oy). How- 
ever, for EMP the form of the equations representing 
the data presented in figures 4-2 through 4-10 is too 
complicated for functional description. Accordingly, 
an alternate method is required, as described in f 
below. 


YIELD = Wy Mt 


Gg = 10-2 mho/m 
OVERPRESSURE = P 
so 
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Figure 4-11. Uncertainty Bands (o = 1) of the B, Fourier Ampli- 
‘tude at the Air/Ground Interface (AF WL, 1972) 
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Figure 4-12. Uncertainty Bands fo = 1) of the Ep Fourter Ampll- 
tude at the Air/Ground Interface (AFWL, 1972) 
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Figure 4-13. . Uncertainty Bands (co = 1) of the E, Fourter Ampli- 
tude at the Air/Ground Interface (AFWL, 1972) 


f Select a mean overpressure, mean weapon yield, 


mean ground conductivity, and mean frequency - 


(DNA, 1972; AFWL, 1972). For example, in figure 4-8 
select yy = Po = Panos ¥9 =We= Wi, ¥3 = Gg = 10-2 
mho/m, and ¥, = f = f,. This defines point A, on 
figure 4-8 in which x = iB, = IBalo- Vary y; within 
the limits, say of P01 and P.,,, which produces a cor- 
responding change in X (represented by the points B 
and C in figure 4-8). Thus, in equation 2-16, for the 
first independent variable y, and forn = 1, . 


_2 
y= (Pyoo) ; 


ae |B,|.)? 


(8g)? = (Ag)? = (By, | - [Bygl)? 


(8yy)? = (yy) 2 = Pog, - Prog)? (4-1) 
The user must select @y, as the uncertainty he 
wishes to assign to the overpressure variable. 
Therefore, for the first component of uncertainty, 
equation 2-16 becomes 


‘4-12 


. (4-2): 


Equation 4-2 is completed by performing the same 
operation for the variables W, Og, and f (if desired) 
and the summation in equation 4-2 is replaced as 
three individual terms corresponding to the uncer- 
tainties in these three variables. 

g. The uncertainty Q, in equation 4-2 is deter- 
mined from the data presented in figures 4-11 
through 4-14 (AFWL, 1972). Here, Q, refers to the 


_ data scatter obtained from Monte Carlo simulations 


performed by the AFWL. . 

h. Note that the Monte Carlo mean values (denoted 
by the symbol e) differ from the nominal mean values 
(denoted by the solid lines). Using the nominal curves 
as a guide, construct new curves through the Monte 
Carlo means and bounds as exemplified in figure 4-11 
for |By(f)| by the dashed lines. These curves may then 
be used to estimate the uncertainties of the Fourier 
spectra at any frequency. ; 

i. From figure 4-11, for example, which cor- 
responds to figure 4-8, the uncertainty 8, is deter- 
mined from equation 2-18 where 


R= 1(At/Ay + Ay/AL)2 
4 | 


n=1 


Here, R is simply the squared averaged ratios of the 
bounds above and below the Monte Carlo mean. Note 
also that the mean value of |B,| (denoted by A, in 
figure 4-11) was obtained from figure 4-8 and subse- 
quently used in equation 4~2, whereas 2, was ob- 
tained from the mean and scatter bands in the Monte 
Carlo simulation. 


|e (f)| AMPLITUDE, V/m/Hz,: LOG SCALE 


-f FREQUENCY, Hz, 
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: , aaa 14. -Fourter Amplitude of E, Waveform of High- Altitude 
hae Se _ Radiated Signal (AFWL, 1972) 


_j. As noted in 6 above, the currents induced in a 
system should be calculated by including the system 
within the EMP enyironment, 1e., by consideri ing the 

system and free” ‘field environment to be coupled. 
When making EMP calculations using computer pro- 
grams, the uncertainties in the field calculations de- 
pend on the uncertainties in the mathematical model 
and the value of its parameters. With these.calcula- 
tions it is recommended that sufficient data be com- 
puted to allow evaluation of equation 2-15. This 
means that a few auxiliary computations should be 
made to allow for determination of the derivatives in 
that equation as outlined in the previous paragraphs. 
The user will be required to ascertain for himself the 
nonrandom and random errors involved in such 
calculations. 


4-3. High-altitude burst. 


a. Any burst occurring at altitudes above 100,000 ft 
(30,480 m) is a high- altitude burst and may produce 
EMP fields that cover ‘geographical areas as large as 
the continental United States. . 

_ 6. The electric field of the EMP is defined by 


=o] 2 vim (4-8) 


where 1, and ty control the’rise and decay times of 
the pulse and E, is the peak amplitude of the pulse 
(DNA EMP Handbook, GE-TEMP, 1971). It is noted 
that the amplitude of the EMP depends on (1) the 
prompt gamma ray yield, which is a function of 
nuclear weapon yield, and (2) the location of the 
system relative to the detonation point. A typical 
pulse is shown in figure 4-15. The corresponding 
Fourier spectrum is presented in figure 4-14. The 
radiated magnetic field B, is proportional to E, (eq. 
4-3) according to the relation 


B= Wo/m?’ | (444) 


where C, is the speed of light in a vacuum (3.0 x 108 
m/s), 

c. The data referenced in GE-TEMPO (197)) should 
be considered mean values of the free- field 
environment. : 

d. Accounting for the uncertainties of high-altitude- 
burst EMP predictions parallels the procedures 
discussed in paragraphs 4-2 d through i for determin- 
ing surface-burst EMP environments, and should be 
calculated accordingly. 
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Figure 4-15. E -Time Waveform of High-Altitude Radiated Signal 
(DNA, 1972) 
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CHAPTER 5 
FIREBALL 


5-1. Gross phenomenology. 


a. In the explosion of nuclear weapons, a large 
amount of energy (kilotons to megatons TNT 
equivalent) is released in a time much less than a 
millionth of a second and in a small] volume extend- 
ing slightly beyond the weapon. The resulting high- 
energy density produces a high-temperature (millions 
of degrees Kelvin, °K), high-pressure, thermally 
radiating plasma core. This core (the fireball) expands 
hydrodynamically because of the high internal 
pressure and radiates because of the high 
temperature. The radiation, primarily in the x-ray 
region of the electro-magnetic spectrum, is absorbed 
and reradiated by the surrounding air. Temperatures 
in the neighborhood of 1,000,000°K persist for a few 
hundred microseconds. As the fireball expands and 
its temperature drops, radiation diffusion becomes 
less important. Hydrodynamic expansion, which pro- 
duces compression heating of the air behind the radia- 
tion front, overtakes the radiation front, developing 
a strong shock that continues to heat air to 
temperatures in excess of 10,000°K for tens of 
milliseconds. 

b. The variables affecting the potential destruc- 
tiveness of the fireball at a ground surface location 
are the weapon yield and slant range. The larger the 
weapon yield, the greater the temperature and the 
engulfment duration will be at a given range. If the 
burst elevation is greater than the fireball radius, 
then engulfment will not occur and the effects of ther- 
mal radiation will control the design. How to deter- 
mine the effects from thermal radiation is presented 
in chapter 6. 

c. Immersion of the facility or its components in the 
fireball, or ingestion of fireball gases by the facility 
can lead to debilitating thermal loads. Knowledge of 
the fireball size and its time/temperature/density 
characteristics are required to determine the design 
strategy—avoidance or resistance—and the degree of 
protection required. 


5-2. Engulfment radius. 


a. The maximum fireball radius is represented by an 
equation of the form 


 Rinax = C0, Wee eft (5-1) 


where W is weapon yield in kilotons (kt); z is HOB 
in ft; and C), n,, and ny are constants. Within this 
radius, facilities will be engulfed by the fireball. The 
temperatures within the fireball are presented below 
in paragraph 5-3. 


b. The mean radius is-calculated from equation 5-1 
using C, = 250 for a contact burst or 220 for a sur- 
face intersecting fireball, n, = 0.35, and ng = 
8.8 x 10-8. 

c. The uncertainty associated with equation 5-1 is 


2 — 02 2 202 202 zg 
Ran Et 20, + ney + 4 Q° (5-2) 


where - is the uncertainty in the form of equation 
5-1, Qc, is the uncertainty of the coefficient C,, etc. 


“Assume that 2,* 0.2 and Qc, = 0.15 (from equation 


2.16 in which n = 2 and R = 1.86). 


5-3. Temperature and density. 


a. The temperature and density within the fireball 
are shown in figures 5-1 and 5-2 as a function of 


‘range for a family of early times. Fireball 


temperature is shown in figure 5-3 as a function of 
range for a series of later times, and in figure 5-4 as 
a function of time for a family of ranges. 

b. The peak fireball temperature vs. range can be 
represented by an equation of the form 


log T, = A + B log (R/W1’3) + C log*(R/W!8) 


+ D log(R/WY3) (5-3) 


T, = Peak temperature in K 
R = Range from GZ in ft 
W = Yield in kt 
A,B,C,D are constants 
c. The mean peak fireball temperature in equation 
5-3 is given by the following coefficients for mean 
peak fireball temperature: 


A B C D 
-479.2 566.4 -235.4 32.36 


d. The uncertainty associated with equation 5-3 is 
Q2 = Oe 5.3 [ + B? log? (R/W"3) + C2 log4 (R/W3) 


+ D? log® nwt) 22 + [p + 2C log (R/WV3) 
0 
+ 3D log? nw) 4 (Q2 + 92/9) (5-4) 


where &, is the uncertainty of the functional form of 
equation 5-3, Qc. is the uncertainty of the coeffi- 
cients, Qp is the uncertainty of the range. Although 

5-1 
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Figure 5-4. Temperature Vs. Time at High Peak Overpressures (adapted from Brode, 1964) 


insufficient data exist to accurately assess the uncer- 
tainty of equation 5-3 at this time, it is recommended 
that 2, = 0.2 and Qqo, = 0.3. 

e. Duration of temperature ts of interest. Since the 
fireball is like a bubble in the atmosphere, it begins 
to rise, and pulls away from the earth’s surface in only 
a few tenths of a second/kt!. Using a very approx- 
imate model for the effect of buoyant fireball rise on 
the temperature history at a scaled ground range of 
265 ft/kt!/3 (80.8 m/kt!/), it appears that the late-time 
temperature environment will be reduced in approx- 
imately the manner indicated by the decreasing tail 


in figure 5-4. Thus a sharp temperature decrease 
begins after 0.4 or 0.5 s/kt!/3 of exposure, and in 1.5 
to 2.0 s/kt! the air temperature has returned to nor- 
ma] ambient. Temperatures at other ranges would be 
similarly reduced at late times by the same effect. At 
the scaled ranges of 205 and 150 ft/kt/? (62.5 and 
45.7 m/kt /3), temperatures increase even more but 
subsequently show a more rapid drop (at times less 
than 0.1 scaled second) because of thermal radiation 
loss, which becomes significant even before the 
fireball has begun to rise. 
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CHAPTER 6 
THERMAL RADIATION 


6-1. Gross phenomenology. 


a. A thermal-radiation pulse is emitted by the ex- 
panding fireball. The thermal radiation received at 
a ground surface of design interest to this volume 
comprises energy in the ultraviolet, the visible, and 
the infrared spectral bands. The radiated thermal 
power reaches a peak in about 1 s/Mt®-5. The pulse 
lasts, for all practical design purposes, about 10 
s/Mt®-5, 

b. The basic variables affecting affecting the poten- 
tial destructiveness of the free-field thermal-radiation 
pulse are the: 

~ Weapon design 

- Weapon yield 

~ Burst elevation 

- Slant range 

~ Transmittance by the intervening atmosphere 

- Reflectance of clouds and of ground surface 
The larger the weapon yield, the greater the height 
of burst, the smaller the slant range, the greater the 
transmittance, or the greater the reflectance, the 
greater will be the potential destructiveness from 
thermal radiation. 

c. The thermal-radiation pulse can cause fires and 
the degradation of the mechanical properties (possibly 
melting) of an exposed facility or facility component. 
The facility designers must either account for these 
effects in their design or preclude them by burying 
(or otherwise protecting) the facility or facility com- 
ponent. Knowledge of two characteristics of the free- 
field thermal-radiation pulse is required to implement 
design measures: Radiant exposure and radiation 
history. 


6-2. Radiant exposure. 


a. The total amount of radiant energy impinging on 
a unit area of a surface that is facing the fireball is 
called the radiant exposure of the surface. The radiant 
exposure is , 


Q = AWFTF/R? cal/cm? (6-1) 
where W = Weapon yield, Mt 
f = Thermal efficiency of the weapon (5 
below) 
T = Transmittance of the intervening at- 


mosphere (c below) 
F = Reflectance of clouds and the ground 
surface (e below) 
R = Slant range, kft 
A = Constant 
Equation 6-1 is applicable for slant ranges greater 
than about 1100kT°4 for an airburst or 145kT®4 for 


a contact surface burst. The radius of the fireball at 
which the airblast shock breaks away and races 
ahead. The radiant exposure at closer points is 
negligible compared to the subsequent heat load,.in 
the enveloping fireball gases, which is treated in 
chapter 5. 

(1) The mean radiant exposure is obtained from 
equation 6-1 where A = 86,000. 

(2) The uncertainty associated with the radiant 
exposure is 


2 02 2 2 24 02 2 2 - 
Sq = B+ RA + My + QE + Sn + QE + 482p (6-2) 


where &, is the uncertainty parameter associated 
with the functional form of Q, Qy, is the uncertainty 
parameter associated with the weapon yield W, and 
so forth. Assume that 22 + 22 = 0.1. The uncertain- 
ties of the transmittance, iiermal efficiency, and 
reflectance are developed in the paragraphs below. 

b. The thermal efficiency f of a weapon is the ratio 
of the thermal energy radiated by the fireball to the 
total yield of the weapon. The thermal efficiency of 
an ordinary nuclear weapon that is detonated in a 
low-altitude burst is 


f=C + Bhiwi8 0 <h < 2.5 kftMt4s (6-3) 


where 
h = Height of burst, kft 
W = Weapon yield, kt 
B and C are constants 
(1) The mean thermal efficiency of an ordinary 
nuclear weapon that is detonated as a low-altitude 
burst is obtained from equation 6-3 where B = 0.078 
and C = 0.21. Assume that the mean thermal effi- 
ciency of an air burst is 0.4, and of a contact burst 
is 0.2. The thermal efficiency of a shallow-buried burst 
is small compared to that for a contact burst. A so- 
called enhanced weapon (enhanced nuclear radiation) 
exhibits lower thermal efficiencies than those cited 
above. 
(2) The uncertainty associated with the thermal 
efficiency of an ordinary nuclear weapon that is 
detonated as a low-altitude burst is 


2p2 
eiiory TS 
f B  f2w2/3 


(0 + 23a] (6-4) 


Assume, 02, = 0.24 for an air burst and = 0.3 for a 
contact burst. 

c. The transmittance T is the ratio of the radiant 
energy reaching a surface facing the fireball to the 
radiant exposure that the surface would receive if the 
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intervening atmosphere were perfectly transparent 
and there were no reflecting surfaces. The transmit- 
tance for ground-surface targets is 


= (1 + aR/V) exp (bR/V) ,h < 1.5 kft (6-5a) 


~T = exp[(o + ph + yh?) rRV], (6-5b) 


1.5 kft <h < 10 kft 


where 
R Slant range, kft 
V = Daytime visibility, kft 
a, b, T, a, B, and y are constants 
Equation 6-5 is based on the assumption that the at- 
mosphere is cloudless, but not necessarily clear of 
haze, fog, smoke, or dust. ° 
(1) The mean transmittance is obtained from 
equation 6-5 where a = 1.9, b = -2.9, T = -84.5, 
a = 3.0, B = -0.286, y = 0.01. Equation 6-5 is plotted 
in figure 6-1. The discontinuity in the transmittance 
shown in the figure is not physical, but is a conse- 
quence of the uncertainty in the two transmittance 
models for the two ranges of burst height in equation 
6-5. 
(2) The uncertainty associated with the transmit- 
tance for the two ranges of burst heights are 


{I 


R2 exp (2bR/V) 
2 — 02 Se 
QF Ser + vey? 
2 2 
(Q5 + 22), 
h < 1.5 kft 


Z 
E + aR/V) + | 


(6-6a) 


Ee) 
to 
| 


= Q2+ e E +: Bh + yh?P (2 + 02) 


+ h%™(B +2yh)? 22, 


15<h<10kft — (6-6b) 


where ¢ is the uncertainty parameter for the func- 
tional form of equation 6-5, Qp is the uncertainty 
parameter for the slant range, etc. Assume 9, = 0.1. 
The uncertainty of the visibility is calculated 
subsequently. 

d. The U.S. Weather Bureau defines day-time 


visibility as the distance at which a prominent dark _ 


object, located near the ground and silhouetted 
against the sky, is visible and recognizable. Night- 
time visibility is defined as the distance an unfocused, 
moderately intense light source is visible. The day- 
time visibility for generic cloudless atmospheric con- 
ditions is given in table 6-1. Assume that the day- 
time and night-time visibilities are identical for iden- 
tical atmospheric conditions. 
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(1) Establish the facility-site, relative-frequency 
visibility function using visibility statistics obtained 
from a local weather station or from personnel at the 
site vicinity. For the seven classes of visibility 
delineated in table 6-1 the facility-site, relative- 
frequency function ¢ is constructed as shown in figure 
6-2 such that.the sum of $95 + $45 + ... $39 = 1.0. 
The subscripts refer to the seven classes of visibility. 
A typical relative-frequency factor is depicted in 
figure 6-2. Convert the weather station’s night-time 
visibilities to equivalent day-time visibilities using 


V = 20 Me kft (6-7) 


ight 

Calculate the mean visibility V as shown in figure 
6-2. Use a base period of one year for calculating 
visibility. 
(2) The uncertainty associated with the visibility 
is ; 
Q2, = Ae poe, (6-8) 
where 

Ay = Coefficient of variation of the mean of 
V relative to its true value 
Coefficient of variation of V, as il- 
lustrated in figure 6-2. 
Assume that Ay = 0.17. . 
' (3) The effect of prior bursts can be important. 
Bursts below about 2900 ft/Mt?-4 create a dust cloud 
that will shroud the facility for many hours. Consider 
the effect of the dust cloud when calculating the 
visibility for each burst of a multiple-burst attack, ex- 


I 


dy 


cept the first burst (chap. 13). 


e. The reflectance F is the ratio of the radiant 
energy reaching a surface facing the fireball when 
there are clouds and/or a reflecting ground surface, 
to the radiant exposure that the surface would receive 
in the absence of clouds and a reflecting ground sur- 
face. The reflectance is 


(6-9) 


PF =F ca F cb F, 
where 
F,, = Reflectance of a layer of clouds above 
the burst 
Fy, = Reflectance of a layer of clouds below 
the burst 
Py = Reflectance of the ground surface 


(1) The mean reflectance is obtained from equa- 
tion 6-9 where the mean values of F,,, F.,, and F, 
are obtained from f g, h below 

(2) The uncertainty associated with the reflec- 
tance is 
#08 _ (6-10) 


ca cb 
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Figure 6-1. Mean Transmittance for a Ground-Surface Target (SAI, 1971) (2 of 5} 
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Figure 6-1. Mean Transmittance for a Ground-Surface Target (SAI, 1971) (3 of 5) 
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Figure 6-1. Mean Transmittance for a Ground-Surface Target (SAI, 1971) (4 of 5) 
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Very clear 


Clear 


Light haze | 


Medium haze 


Heavy haze 


Thin fog 


Light fog 


Table 6-1. Day-time Visibility for Generic Cloudlessness (DNA, 1972) 


Description 


This condition rare except at 
high-altitude locations. 


Sky deep blue. Shadows distinct, 
dark. 


Sky white; dazzling near sun. 


Shadows visible, gray. 


Sky bright grayish-white. View 
sun without serious discomfort. 
Shadows visible but faint. 


Sky dull gray-white. Sun's disc 
just ‘visible. Shadows barely 
discernible. 


Sky light gray with maximum 
luminance around sun. Sun's disc 
not visible, no shadows. 


Sky dull gray with maximum 
luminance at zenith. 
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f When the burst occurs below a cloud layer, the 
reflectance iS 


Fig = A+ Bog (6-11) 


where A and B are constants and $,, is the fraction 
of a year that clouds appear in the site vicinity at a 
height above burst elevations of interest. 
(1) The mean clouds-above reflectance is obtained 
from equation 6-11 where A = 1.0 and B = 0.5. 
(2) The uncertainty associated with the clouds- 
above reflectance is 


Qe = Az 0:28 tall bu) (6-12) 
ge A OB bes 

where AF oa is the coefficient of variation of F., 

relative to its true value. Assume that Ap, = 0.1. 
g. When a layer of haze or clouds is below the burst, 

but high enough not to affect visibility, the reflectance 

F.,, is as shown in table 6-2. For this condition, use 

VY = 16 mi in the calculation of transmittance. 

(1) Calculate the mean value of the clouds-below 
reflectance as illustrated in figure 6-3. Use a base 
period of one year for calculating the mean value of 
Fy- The method of calculating the mean reflectance 
parallels the procedure described in d above. 


Table 6-2. Reflectance F, for Layer of Clouds or Haze Below the Burst (DNA, (1972) 


| tyre | Description Reflectance 


Very clear 


Clear Sky deep blue. 


dark. 


Light haze 


This condition rare except at high- 
altitude locations. 


Shadows distinct, 


Sky white; dazzling near sun. 


Shadows visible, gray. 


Medium haze 


Sky bright grayish-white. 


View 


sun without serious discomfort. 
Shadows visible but faint. 


Heavy haze 
just visible. 
discernible. 


Light cloud 


Sky dull gray-white. 
Shadows barely 


Sun's disc 


Sky light gray with maximum 
luminance around sun. 


Sun's disc 


not visible, no shadows. 


Medium cloud 


Sky dull gray with maximum 


luminance at zenith. 


Heavy cloud 


Sky dark gray; brightness pattern 


gives no indication of sun's 


position. 


Dense cloud 


The low luminance level suggests 


the approach of nightfall. 


Murky. 
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(2) The uncertainty associated with the clouds- 
below reflectance is 


Q2 =A2 + & 6-13) 
Fob Fob Fob 
where 
4p = Coefficient of variation of the mean 
value of F.,, relative to its true value 
dp, = Coefficient of variation of F,,, as 


calculated in figure 6-3 
Assume that Ap, = 0.40 - 0.35 Fy, 

h. The ground surface pelleclence 3 in n the facility- site 
vicinity is given by 
Fo= Ky Koo, (6-14) 
where $5 is the fraction of the year that the ground 
surface appears bright..Snow and light- colored soil 
qualify as bright surfaces. 

(1) The mean ground-surface vanlebianes is 
obtained from equation 6-14 where K,; = 1 and 
Ky = 0.5. 

(2) The uncertainty associated with ground- 
surface reflectance is 


A Az + 0.25 4g - 4) (6-1) 
g gE (14+05 4.) 

where Ay,, is the coefficient of variation of F, relative 

to its true value. Assume that Ap, = 0.1. 

(3) Dust created by prior bursts will settle out and 
decrease the reflectance of a bright surface. Consider 
this effect when calculating the reflectance for each 
burst of a multiple-burst attack, except the first burst 
(chap. 13). 


6-3. Radiation history. 


a. The history of the radiation energy is relatively 
complex, exhibiting several peaks and valleys. The 
details of the history are strongly influenced by the 
weapon yield, and the burst elevation. For purposes 
of design, the simple fractional radiant-exposure 
history shown in figure 6-4 can be used in place of 
the radiation history. The functional form of a curve 


fitted to the fractional radiant-exposure history of 
figure 6-4 is 


Q.=Q (1 - exp [atit, + p (tt od] 


cal/em?, 1 < tit, < 10 (6-15) 
where 

Q = Total radiant exposure, cal/em? 

t = Time, seconds 

t, = Normalizing time, seconds (b below) 


a, §f are constants 
(1) The mean fractional radiant-exposure is 


obtained from equation 6-16-where a = -0.5 and 


B = 0.33. 
(2) The uncertainty Asaoriated with the fractional 
radiant-exposure is 


2 - 924 Q2 _Q) 2]2 
QB = 92408 +0-9) [atit, + 2p) 


(6-16) 
Assume Q, = 0.1. 
b. The normalizing time t, is 
t, = AWe{1 + Bh + Ch?]? sec (6-17) 
where i‘ 
W = Weapon yield in Mt 
h = Height-of burst, kft 


A, B, C, ¢, o are constants 
(1) The mean.normalizing time is obtained from 
equation 6-18 where A = 0.84, ¢ = 0.431, B = -0.0162 
C = 7.2 x 105 and o = 0.764. The mean normalizing 
time can also be obtained from figure 6-5. 
(2) The uncertainty associated with the normaliz- 
ing time is 


2 —~ o2 4 2 02 2 
oF =Qit+e of, + othe 


i) 


B + 2Ch i a 
1+Bh+cCh2] > 
. (6-19) 


where Qy is the uncertainty parameter for the yield 


‘and Q,, is the uncertainty parameter for the HOB. 
Assume &, 


= 0.17. 


Q./Q FRACTIONAL RADIANT EXPOSURE 


k 6 
t/t, NORMALIZED TIME 


Figure 6-4. Fractional Radiant Exposure 
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CHAPTER 7 
AIRBLAST 


7-1. Gross phenomenology. 


a. Hydrodynamic expansion of bomb debris and 
fireball produces an expanding shock wave in the sur- 
rounding air. Interaction with the ground surface af- 
fects the detailed structure and geometry of the ex- 
panding shock wave. 

b. A contact burst at a scaled elevation no greater 
than 5 ft/kt}/3 (1.5 m/kt!/3) produces a hemispherical 
shock front expanding supersonically. At scaled 
heights of burst between 5 and 160 ft/kt1’3 (48.8 
nV/kt’3), the incident spherical shock reflects from the 
ground (fig. 7-1); the reflected signal, traveling at a 
higher speed in air heated by the incident shock, even- 
tually overtakes the incident shock, forming a Mach 
stem that connects the ground with the.triple point 
(which is the intersection of incident, reflected, and 
Mach stem shocks). The dynamic movement of the air 
behind the Mach stem causes a second shock to form, 
producing double shock phenomena (Brode-Lewis, 
1975). Thermal-radiant heating of the ground ahead 


INCIDENT 


SHOCK ——~__ 


of the Mach stem causes the formation of a precursor 
shock, which trails back and merges with the Mach 
stem. The aspects of airblast important to the 
designer are the overpressure and dynamic pressure 
pulses. These are often described by the principal 
parameters: peak pressure, positive phase duration, 
and positive phase impulse (the time integral of the 
pressure history). Overpressure and dynamic pressure 
histories are also required; they are summarized in 
paragraphs 7-2 and 7-3. 

c. The basic variables governing airblast are: 

- Weapon design 

- Weapon yield 

~ Range from detonation 

~ Burst elevation 

- Surface conditions 

- Surface contour 

- Structural configuration 
The larger the weapon yield, the greater the height 
of burst up to the limit of the Mach, stem, or the 


"Rip, 
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SAS fe SECOND SHOCK 
MACH STEM _<. see 
LY ~~ 
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Figure 7-1. Composition of Shocks Due to Nuclear Airburst (Brode, 1964) 
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shorter the range, the greater will be the potential 
destructiveness of the airblast. The influence of the 
other variables is more complex; changes in weapon 
design and terrain conditions may amplify or at- 
tenuate the local airblast environment. 

d. At the ground surface or in exposed configura- 
tions, structures can be dealt a crippling or fatal blow 
by the airblast. The facility designer must either ac- 
count for the effects in the design or prevent the struc- 
tures’ direct exposure to airblast. To determine the 
degree of resistance necessary and to support selec- 
tion of the design approach, knowledge of the airblast 
waveform, peak pressure, positive phase duration, 
and positive phase impulse are required. 


7-2. Free-field peak overpressure. 


a. Surface to moderately high airburst. The peak over- 
pressure P., depends on scaled HOB ‘height of burst), 
scaled range from GZ (ground zero), and terrain 
characteristics. In general, the peak overpressure P,,, 
may be defined as a piecewise function of scaled range 
for a given scaled HOB, as indicated in the following 
equation: 


= R\" | 
P,, = Cy | wi psi, Ry. 1 <R<R, 


p=1,2,...,N (7-1) 
where R is the range from GZ in ft and W is yield in 
kt. The coefficients C, and exponents n, depend on 
scaled range, scaled HOB and the terrain characteris- 
tics. 

(1) Ground surfaces are characterized as therm- 
ally “near-ideal” (in which a precursor shock is not 
likely to form) or thermally ‘‘nonideal” (in which a 
precursor shock is likely to form). Table 7-1 presents 
a compilation of ideal and nonideal surfaces. These 
categories are used to define the parameters of equa- 
tion 7-1. 

(2) The mean values to be used in equation 7-1 
are presented in table 7~2 for ideal and nonideal sur- 
faces. The corresponding data from which the table 
was prepared are presented in figures 7-2 and 7-3. 
These data are also presented in figures 7-4 through 
7-10, which should be used directly for selecting 
mean values. 


Table 7-1. Thermally Near-Ideal and Thermatly Nonideal Surfaces (DNA, 1972) 


Thermally Near-Ideal 
{precursor unlikely) 
Water 


Ground covered by white smoke 
layer 


Heat-reflecting concrete 
Frozen tundra 


Ice 


Packed snow 


Soil with sparse vegetation 


Commercial and ‘industrial areas 


Thermally Nonideal 
(precursor may occur) 

Desert sand 

Coral 

Asphalt 

Surfaces with low, thick vegetatio. 


Surfaces covered by a dark smoke 
layer 


Dark colored rock 


Most agricultural areas; 
residential areas in cities 


Dry soil with sparse vegetation 


*Note: 


Scaled HOB 


-ftskel/3 


Table 7-2, Parameters for Peak Overpressure 


Scaled Range 


ensue 2/5 


330 


1000 to 


270 
330 
1220 


to 


ie) 


290 
500 
1350 


150 
340 
1000 
240 
750 
1200 


210 
480 
1650 


140 
360 
880 


300 
600 
1500 


300 
“800 
1920 


550 
2500 


170 to 
to 


> 6400 


to 5300 


330 
1000 
3800 


330 
1220 
4400 


500 
1350 
4900 
$20 


1520 (158 to 463) 
5600 


1500 
6000 


1920 


2500 


Compute values from Figures 7-2 or 7-3. 
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Figure 7-2. Peak Overpressure at Thermally Near-Ideal Surfaces 
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Figure 7-3. Peak Overpressure at Thermally Nonideal Surfaces 
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Figure 7-7. Peak Ouverpressures at a Near-Ideal or Thermally Nonideal Surface, Very Low Overpressure Region (DNA, 1972) 
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Figure 7-10. Peak Overpressure at the Surface from a Contact Surface Burst (DNA, 1972) 


(3) The uncertainty associated with the peak 
overpressure Is: 


2 


n 
92 = 024 02 292 4 “Poe 
ms SF + =C, + ni QR + 7 Qo, (7-2) 


where 9, is the uncertainty associated with the func- 
tional form of p,,, 2c, is the uncertainty of the coef- 
ficient, and so forth. For practical use, the major 
uncertainty in equation 7-1 as a representation of ex- 
‘perimental and theoretical data is in the coefficient 
C,,. It is recommended that 2 = 0.05 be used and that 
Qc. = 0.2 (which assumes that the data are log- 
normally distributed and that the bounds of +2 stan- 
dard deviations are located between 3/2 and 2/3 times 
the mean value). 

6. High-altitude bursts. Peak overpressures from 
high-altitude bursts (large HOB) are obtained from 
figure 7-11, which shows peak pressure as a function 
of scaled slant range. These data may be used at 


7-12 


points other than at surface zero if the target remains 
in the region of regular reflection, that is, if the burst 
point is no more than 30 deg (0.52 rad) off the vertical. 

(1) The equation of the data presented in figure 
7-11 (construed to be the mean value) can be 
represented by: - 


P_ = 104 + Blog (R/W?}") + C logtR/W2) + D log(R/W2°) 
sO . : 
(7-3) 


where P,, is the peak overpressure in psi, R, is slant 
range in ft, W, is the effective yield in kt, and A, B, 
C, and D are constants. The effective yield is given 
by the Blast Efficiency Factor; 


Ww. 
—£=¢@+ Bh + yh? (7-4) 
W B Y 


where W is the actual yield in kt, h is the HOB in 
kft, and a, #, and y are constants. Equation 7-4 is 
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Figure 7-11. High-Altitude Bursts, Peak Reflected Overpressure at Ground Zero as a Function of Scaled Slant Range (DNA, 1972) 


plotted in figure 7-12. It is preferable to use data — 


directly from figure 7-11 and 7-12 rather than from 
equations 7~3 and 7-4 to calculate the mean over- 
pressure and the mean effective yield. 

(2) The uncertainty of the overpressure at the 
ground surface due to a high-altitude burst as 
represented by equation 7-3 is: 


OF Bt 5.3[A2 + B2logR/W2) 


+ ClogtR/W13) + DlogS(R/W2) ] 9? 
+ [B + 2C log(R,/W)’3) 
Q2 
+ BD log’R,/W23)]” (Q2 + —2) (7-5) 
- s 9 


where @, is the uncertainty in the form of equation 
7-3, Sp, is the uncertainty of the slant range, etc. 
Assume that 2, = 0.1 te account fer inaccuracies in 
fitting the function in equation 7-3 to the curve in 
figure 7-11. Also let 


1 , 
p——p'[e2 a p*h2 Be Ph ]22 0, 
a+ Bht+ re 


+ n>[ 5 + yh}? ge (7-6) 


The values of the constants in equations 7-3 through 
7-6 are A = 12.0, B = -6.08, C = 0.976, and 
D = ~0.0608; for h < 40 kft, a = 1.0, 8B = 0.0015, and 
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y = 0; for h > 40 kft, a = 0.956, B = 0.00101, and 
y =-3.55 x 10-5. The uncertainty of Qcor is ob- 
tained by applying equation 2-19 to the data in figure 
7-12 where Ly and L, are the upper and lower 
bounds of the Blast Efficiency Factor at any particular 
HOB. Consider that 2c, = 0.13 to account for the 
uncertainty in the data comprising figure. 7-11 
(derived from equation 2-18 in which VY R = n = 2). 


7-3. Free-field peak dynamic pressure. 
a. At GZ the ground-surface dynamic pressure is zero 


for all HOB. Outward from GZ, the peak dynamic . 


pressure rises and then falls with increasing range. 
Unlike overpressure, the dynamic pressure is 
relatively insensitive to the surface properties and 
depends more on the dust loading. Peak dynamic : 
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pressure, a complicated function of HOB and range, 
is not easily represented by simple linear relation- 
ships, regardless of the manner in which the data are 
presented. Figures 7-13 through 7-15 present 
dynamic pressure curves in their simplest forms. In 
general, these curves may be Fepresented by 


R \2 R 
log Q =A, (2) + B. (| + C,, 
R,,<R<R, i=1,2,...,N (7-7) 


where Q is the dynamic pressure in psi, R is the range 
from GZ in ft, W is yield in kt, and A, B, and C are 
coefficients to be defined in c below for various HOBs 
within certain scaled ranges. 
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Figure 7-13. Dynamic Pressure at a Perfect Surface (DNA, 1972) 


7215." 


T™<\ 5-858-2 


SCALED RANGE, m/ktl/ 
0 100 200 300 400 


SCALED HOB = 400 ee (121.9 m/kt!/3) 


f a (121.9) 


SURFACE 


Q, DYNAMIC PRESSURE, psi 


"G00 (182.9) 


0 200 400 600 800 1000 1200 1400 


r/w!/3, SCALED RANGE, Ft/kt'/3 . 
Figure 7-14. Dynamic Pressure at a Light Dust Surface (DNA, 1972) , 
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Figure 7-15. Dynamic Pressure at a Heavy Dust Surface (DNA, 1972) 


b. It is recommended that the mean values for 
dynamic pressure be taken directly from figures 7-13 
through 7-15 or from the source data presented in 
figures 7-16 through 7-18. 

c. The uncertainty associated with equation 7-7 is 


Q2 = 


2 2, 1/3)4 2 1/3)2 2! O2 
2 = 92+ 5.3)[A%R/W!S)4 + BAR/WI9? + C2] 2 


Q2 
ah 
(7-8) 


where &- is the uncertainty in the form of equation 
7-7, Qo, is the uncertainty of the data in figures 7-13 


+ (RW)? [2aR/w"’) + BP [92 + 


through 7-15, 2y, is the uncertainty of the yield. Due 
to the complexity of the curves in figures 7-13 
through 7-15, the coefficients A, B, and C have not 
been calculated in advance. They can be determined 
for any range of the data by fitting equation 7-7 to 
the data within the range of interest. For example, 
if the curve corresponding to the scaled HOB of 400 
ft/kt1/? in figure 7-14 is fitted by equation 7-7 in the 
scaled range 570 < R/W!8 < 1200 ft/kt1/3, the approx- 
imate regression equation will be specified when 
A = 1.021 x 10-®, B = -3.485 x 10-3, and C = 4.013. 
Credible estimates of the values of 2, are not cur- 
rently available from the library of data. However, 
for design purposes assume 2; = 0.1 and 2, = 0.3. 
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Figure 7-16. Peak Dynamic Pressure at an' Ideal Surface (DNA, 1972) 
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7-4, Overpressure positive phase duration. 


a. The positive phase duration, that is, the time that 
the overpressure remains positive, is relatively in- 
sensitive to HOB or range. Reasonably accurate 
mathematical representations of the positive phase 
durations presented in figures 7-19 and 7-20 may be 
achieved by interchanging the variables on the 
ordinate with the variables on the family of curves 
and replotting the data as shown in figures 7-21 and 
7-22. Equations of the latter forms of the positive 
phase durations are 


+ 


D5 
— =C, _ + Cp, log — 


ma Rj <R<R, 


ws 


ee ere | (7-9) 


for a scaled HOB of 200 ft/kt!/3 (61 m/kt!’8) or 


+ 


Riek <P: 


oP = Ci a 
wi3 oj Di wi 


i=1,2,...,N (7-10) 
for other selected scaled HOBs. In equations 7-9 and 
7-10, D~ is the positive overpressure duration in sec, 
W is yield in kt, R is range in ft, and C, and Cp are 
constants. 

b. It is vevommiended that the mean values of the 
positive phase durations be extracted from figures 
7-19 and 7-20 by interpolating the plotted data 
presented in table 7-3. 

c. The uncertainty of the positive phase duration 
for equation 7-9 is to be computed from: 


Q2 = Q2+ we C2 + C2 log? R/W3! Q2 
D+ “f D+2 fo) D 8 c 
p + 

p 
202 
+ 0.189 CEee 


2 
+ dey (0.494 — tog RAW) -c,]" 02 


(7-11) 
and for Equation 7-10: 
‘ 
2 — O2 a 2/802 20:2 2 
Q.= Qe + -[w c2 + re2, Jos 
p 
P 
W2 202 202 02 
+ —— C02, +R Cy Qe (7-12) 


TM 5-858-2 


where 2; is the uncertainty. of the functional form of 
equations 7-9 and 7-10, 2, is the uncertainty of the 
coefficients, and 2p is the uncertainty of the scaled 
range. The value of the constants in equations 7-9 
through 7-11 are presented in table 7-3. Assume 
Q-= 0.2 and Q, = 0.1 for D3/W'3 < 0.25 s/kt! or 
2,= 0.17 for “Dewi > 0. 25 s/kt’/3 (which is 
predicated on the use of equation 2-17 where n = 2 


and R = 1.4 or R = 2, respectively). 


7-5, Dynamic positive phase 


duration, 


pressure 


a. Dynamic pressure positive phase durations are not 
as readily defined as the over-pressure positive phase 
durations. The usual practice is to modify the positive 
phase of the overpressure by certain factors. Thus, the 
dynamic pressure positive phase is approximated by 


Dt 
wis 


RAW, Ry 


<R<R, 


(7-13) 


from figure 7-23, which is recast from figures 7-24 
and 7-25 where C, and Cp are constants. The 
dynamic positive phase duration Djhas units of sec, 


'W is the yield in kt, and R is the range in ft. 


b. The mean values of the positive phase should be 
extracted from figures 7-24 and 7-25 by interpolating 
the plotted data. 

c. The uncertainty of the dynamic pressure positive 
phase duration is computed from 


1 
2 =~ 02 2/3 72 2 72] 02 
Qe, Qe + =f C* +R ce Q 


WEE Roe a ep Ales 
+—— CQ +R Ch QF (7-14) 


where 2 is the uncertainty of the functional forms 
of equation 7-13, Q, is the uncertainty of the coeffi- 
cients of the equation, Qp is the uncertainty of the 
range, etc. The values of the coefficients are presented 
pete seen eto ad Oro 

< 0.35 s/kt!’3, and Q, = 0.17 for Ds/wis > 0.35 s/fktU3 
(from equation 2- 17 where n = 2 and R = 1.3 and 
R = 1.5, respectively). 
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Figure 7-21. Duration of Positive Overpressure at the Surface 
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Figure 7-22. Birwion of Positive Overpressure at the Surface 
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Table 7-3. Parameters for Positive Phase Duration 


Tels (m/kt !/ ft/kt i! (m/kee? 


300 to 750 (91 to 
750 to 2150 (229 to 


2.50 x 10% 
1.48 x 10° 


440 to 860 (134 to 
860 to 2450 (262 to 


3.44 
3.79 


B39: 107 
1.37 x 107" 


450 to 850 (137 to 
850 to 2075 (259 t 


350 to 780 (107 to 
780 to 1650 (238 to 


1.49 


4.20 
2.41 
4.03 


380 to 500 (116 to 
' 500 to 1300 (152 to 
1300 to 2300 (396 to 


3.00 
742: 


150 to 600 (46 to 
600 to 2150 (183 to 


U.S. Army Corps of Engineers 
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Table 7-4. Parameters for Dynamic Pressure Positive Phase Duration 


Petsetl/S — anserl/3) ecsutll? m/e!’ 5 @ oF 


Near-Ideal Surface 


350 to 550 (107 to 168)" 
550 to 1100 (168 to 335) 
1100 to 2150 (335 to 655) 


450 to 600 (137 to 183) 
600 to 1100 (183 to 335) 
1100 to 2450 (335 to 747) 


0.000040 
. 00010 


450 to 650 (137 to 198) 
650 to 1250 .(198 to 381) ° 
1250 to 2000 (381 to 610) 
2000 to 3100 (610 to 945) 


. 00040 
. 000033 
.00011 
.000055 


Nonideal Surface 


0.00053 
0.000042 
0.00012 


. 00080 
0.000029 
. 000067 


0.00053 
000054 


U.S:- Army ‘Corps of Engineers 


7-6. Pressure waveforms. - Effect of rain, fog, or snow : 
a. General. Under ideal conditions, i.e., airblast pro- - Terrain features 
pagating across snow, ice, or water, the overpressure - Weapon design 
and dynamic waveforms can be characterized by the This document addresses only the precursor forma- 
curves presented in figures 7-26 and 7-27, where the tions since this probably will predominate for surface 
peak pressure and positive phase durations are ob- structures. Most structures will be located below the 
tained from the data presented in paragraphs 7-2 triple point and will be designed for worst-case con- 
through 7-5. Under nonideal conditions, such idea- ditions (clear day) over the actual terrain. The terrain 
lized waveforms will not apply, since a number of fac- (topography) must be considered in determining the 
tors act to distort and modify the blast shape. Among actual waveforms. Exposed structures locally convert 
these are: free-field dynamic pressure to overpressure; the A 
~ Location of the precursor formation region degree and areal extent of conversion depend on size 
- Height of a structure relative to the location of and orientation of the structures (DNA, 1972; GE- 
the regular and Mach reflection regions TEMPO, 1972 reports.) 
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Figure 7-26. Positive Overpressure Waveforms for an Ideal Shock Wave in a Standard Sea Level Atmosphere (DNA, 1972) 
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Figure 7-27, Positive Dynamic Pressure Waveforms for an Ideal Shock Wave in a Standard Sea Level Atmosphere (DNA, 1972) 
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6. Precursor effect—overpressure. Waveforms are 
classified in terms of the state of precursor formation 
as the five general types shown in figures 7-28. The 
type of waveform for a given situation is to be 
estimated from the data presented in figure 7-29. The 
peak values of the pulse, the positive phase duration 
and the impulse should be obtained from the data 
presented in paragraphs 7-2, 7-4, and 7-7 (GE- 
TEMPO, 1972). 

c. Precursor effect—dynamic. pressure. Dynamic 
pressure waveforms are also categorized as the five 
general types as shown in figure 7-30. The ap- 
propriate type for a given problem is estimated from 
figure 7-29. Actual peak dynamic pressures should 
be obtained from paragraph 7-3, and the positive 
phase should be obtained from sperAreDe 7-5 (GE- 
TEMPO, 1972). 


7-7. Overpressure ‘agus: 


a. Overpressure impulse at the surface is primarily 
a function of.slant range. Reasonably accurate equa- 
tions of the impulse. data presented in figures 7-31 
and 7-32 (which show parametric representations of 
impulse-as functions of HOB and ground range) can 
be represented by 


Z R R 


wi3 wis 


R 
+ Cy exp (Cs ae (7-15) 


where I is the impulse in psi-s, W is yield in kt, and 
R is range in ft. The symbols C, through C, denote 
constants. 

b. The mean values for overpressure impulse can 
be obtained directly from figures 7-31 and 7-32 for 


TM 5-858-2 


near-ideal and nonideal surfaces, respectively, or from 
figures 7-33 and 7-34, or from equation 7-15 in con- 
junction with the data in table 7-5. It is recommended 
that the mean values be obtained from figures 7-31 
and 7-32 by interpolation. 

c. The uncertainty of the overpressure impulse cor- 
responding to equation 7-15 is 


Q2= O24 


1 
Ez \cew29 + CoRA+ C2W29 exp(2C5R/W3) 


+ Cowes exp(2C,R/W!8) Q2 


ws | c, _ Cask exp(O;R/W a) 
9 a Wwe 


+ Cg exp(C,R/W!8) 


C4CgR exp(CgR/W12) 


2 
13 2 
= + Cy exp(CgR/W | Qe 


. 2 


where ®, is the uncertainty of equation 7-15, Q, is 
the uncertainty of the coefficients, 2p is the uncer- 
tainty of the eroune range, and so forth. Assume 
Qe= 0 and Q,= 0.1 (using equation 2-17 where 
n= 2 and R = 1.5). 


7-33 


TM 5-858-2 


TYPE | 


PREPRECURSOR 


TYPE |} 


STAGES OF DEVELOPMENT OF PRECURSOR . : 
TYPE HII. | ae | 
STAGES OF EVOLUTION OF FULLY DEVELOPED PRECURSOR 


TYPE IV = ee _ c —— — 


STAGES OF DECLINE OF PRECURSOR 


POSTPRECURSOR 


Figure 7-28. Typical Sequence of Overpressure Waveforms over Thermally Nonideal Surfaces (DNA, 1972) 
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Figure 7-30. Typical Sequence of Dynamic Pressure Waveforms over Thermaily Nonideal Surfaces (DNA, 1972) 
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. Figure 7-33. Positive Overpressure Impulse at a Near-Ideal Surface different burst heights 
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Figure 7-34. Positive Querpressure Impulse at a Nonideal Surface different burst heights 
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CHAPTER 8 
AIRBLAST-INDUCED GROUND SHOCK 


8-1. Gross phenomenology. 


a. At near-surface locations, airblast-induced ground 
shock can cause outright structural failure of facilities 
or the failure or malfunction of internal systems. To 
achieve survival, whether by burial or other harden- 
ing techniques, the facility designer must have 
knowledge of the intensity of the ground-shock en- 
vironment and of the waveforms. , 

b. The basic variables affecting the intensity of the 
ground-shock environment are: 

- Weapon design 

~ Weapon yield 

- Burst elevation 

- Slant range 

- Earth material properties 
The larger the yield, the lower the HOB, the smaller 
the ground impedance, and the closer the facility to 
the surface, the more severe will be the airblast- 
induced ground-shock environment. 

c. The airblast signal produced by a near-surface 
nuclear explosion causes ground-shock signals to be 
formed as the shock front expands. Near the surface 
in the high overpressure region, the ground-shock 
signal is a sequence of events represented as the com- 
bination of (1) a downward, outward motion from in- 
itial loading by the overhead airblast overpressure, 
and (2) a later upward, outward motion from release 
of the overhead overpressure and from motion in- 
duced by overpressure closer in but delayed by trans- 
mission through the ground. The first resembles the 
airblast overpressure in pulse shape and duration, the 
pulse shape modifying with depth. The second con- 
tribution, termed a “run-up” contribution because the 
effects are delayed by travel time through the 
medium, comes from the intense overpressure occur- 
ring in the region near GZ and the subsequent com- 
plex of motions known as directly induced ground 
shock. - 

d. At more distance points, outrunning motion 
occurs in which the signal in the ground outraces the 
air shock on the surface. The motion is characterized 
by an upward and outward surge and often exhibits 
an oscillatory nature. In addition to these contribu- 


tions, there may also be crater-induced motion, which 


is most pronounced near the crater and diminishes 
with increasing range. 

e. In the region between the ground surface and a 
depression angle of about 30 deg below the surface, 
the response is arbitrarily termed airblast-induced 
ground shock regardless of the origin of the motion. 
See figure 9-1 to visualize delineations made in thi 
volume. 


8-2. Ground Response. 


For convenience, airblast-induced ground response is 
usually segregated into two regimes—superseismic 
and outrunning. In the superseismic region, the speed 
of the overpressure shock front is greater than the 
shock speed in the ground. The ground-shock signal 
appears to be attached to and dragged along by the 
overpressure shock front traveling across the surface. 
In the outrunning region, the ground-shock signals 
arrive first, “outrunning”’ the overpressure shock 
front. 


8-3. Superseismic peak soil response. 


a. Peak soil acceleration. Peak horizontal and vertical 
accelerations are directly related to the peak over- 
pressure at the surface, modified by the soil proper- 
ties and the depth. The vertical component is given 
by Sauer (1964) and Crawford et al. (1974) as 


80 Zn 


(8-1a) 


ay = k, 
oY 


or 


k, Zn 1/3 1/3 
= 07 419A + B log RW"? + C log? RW 


ay = 
oY 
where 
ay = Vertical soil acceleration, g 
P,, = Overpressure, psi 
C, = Soil seismic sound speed, ft/s 
y = In situ soil density, lb/ft? 
Z = Depth, ft 
R = Range, ft 


W> = Yield, kt 

A, B, C, n, k, = Constants 
The horizontal component of soil acceleration may be 
taken as 


Cc 
Vue 5 C,< UY 2 


ay = 


A ay . CAS UY 2 (8—1b) 
where U = overhead airblast shock speed, ft/s (see fig. 


29 in Brode, 1964) 

b. Peak soil velocity. Peak vertical and horizontal 
velocities are represented by 
uy = 0.11 (Z,W,C,) P9896 (8-2a) 


or 
8-1 
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uy = 0.11 f(Z,W,C,) 1QA' + Blog RWS + C’ log2R/wi - 


and © 
uy = D(R/W18) (8-2b) 
where 
uy = Vertical soil particle velocity, ft/s 
UH = Horizontal soil particle velocity, 
ft/s _ 
f(Z,W,Cy) = ee 
eR 415 Z + 11 C, w¥3 
CL = Compressional wave speed of 


soil under load, ft/s 
_ A’, B’, C’, D, m are constants 

Equations 8-2 are adaptations of the work presented 
by Sauer (1964) and Murrell (1974). Note that the 
horizontal particle velocity does not attenuate with 
depth. ; 

c. Peak soil displacement. The peak vertical and 
horizontal. displacements are (Sauer, 1964) 


v3 P.. - 40 
dy = Z,W,C,) 413 Pot + ogo Swi, in 
; : ° Ba) 
dy = dy/3, in. (8-3b) 


where the variables have been previously defined. 
Equations 8-3 may be converted to range via the 
relationship 

log P,, = A + B log R/W13 + C log? R/W1 (8-4) 


d. Peak soil stress. The peak stress in the soil is 


computed from equations 8-2 via the relationships. 


C 
oy = EY psi (8-5a) 
144 
Oy = —— ey, psi (8-5b) 
-v 
where 
_ ¥ lb sec? 
g .ft# 


g= 32.2 ft/sec? 
v = Poisson’s ratio 
e. Constants in mean and uncertainty calculations 
of peak soil response. For peak accelerations (eqs. 8-1, 
peak velocities (eqs. 8-2), peak displacements (eqs. 


8-2 


8-3 and 8-4), and peak stress (eqs. 8-5), use the 
following constants to calculate mean values: 


1.72 x 10° -04 , 1<Z< 20ft 
k, = n= 

4.36 x 106 -148 , Z > 20 ft 
A = 11.87 A’ = 10.64 
B = -5.26 B’ = -4.71 
C = 0.543 Cc’ 0.487 
D = 2.72 x 105 
m = -2.09 


In calculating uncertainties at fg, h, and i below, find 
Qp,, from equation 7-2 and assume the following 
value: 

0.13 


Q, = 

2, = 0.09 
Q, = 0.24 
2), = 0.138 
Q, = 0.10 


m 
f Uncertainty of peak soil accelerations. Using the 
values from e above, calculate uncertainties of ver- 
tical and horizontal accelerations: 
2 2 2 2 202 
Se + i + ae. or + n*Qr 


or: 


Q2 = 


24 92 24 p22 
ay PEt QE + 22 + n2Qe 


ce) 


+ 


2fQ2 
B + 2C log Rw e + 25 | (8-6a) 


1 
Q2 + Q2 ———— 
PS (U2 — C22 


+ Q24 n2Q2 + 
sO yo. Z 


‘T1402 402 : 
lu Q2, + cick | Cy SUAS 2 


9? or 
aH 


If 


1 
(U2 - 0%? 


24 024 n202 402 2 
Q2 + oe + n?Q2 + U4Q2 + coe | 


| ues 
~ [B + 2C log Rw} [=e + 24) ; 


C,<UN 2 
Q2 Co> UN 2 


(8-6b) ~ | 


g. Uncertainty of peak soil velocities, Using values 
from e above, calculate uncertainties of vertical and 
horizontal peak velocities: 


' 72602 
| ee je +92 +92 
9 CL Z 


& 11115Z+ 110, Ww 
+ 0.80322 
Ps 
Tr 
2 
Q2 = 22+] atta | cE +f 
uy 6 11115Z+11C,w'3) bc. 2 


2 
+ IB + 2C’ log nw:s| 92 


1115 Z aa 


+ 
1115Z +110, Ww 


2 
: 2C log rwis| ore (8-7a) 
2 2 i 2 Qi 
or = 02 + m7Q2 + ; | (8-7b) 


h. Uncertainty of peak soil displacements. Using 
values from e above, calculate uncertainties of ver- 
tical and horizontal peak displacements: 


1115 Z ; 
op -[__s2_ Flog a 
dy kK 11115 Z+110,wis] | cL. 4 
2 
[p280 Z + 11 CLWIS? QE 
1115Z2+11C,Wi3l 9 


5.2P04C,+280P., 2 


Q2 . (8-8a) 
13 P04 C, + 280P,, - ral Poo 
Q? = 22 (8-8b) 
H dy 
These equations may be expressed in terms of range 
and yield by utilizing equation 7-2. 
i. Uncertainty of peak soil stress. Using values in 
e above, calculate uncertainties of vertical and 
horizontal] peak stresses: 


(8-Qa) 


Ore OF + Qh + + 92 
oy uy 

G2" 2. die ' (8-9b) 
OW oy v 


8-4. Superseismic peak rock response. 


a. Peak rock acceleration. The peak horiontal and ver- 
tical accelerations in rock are described by equations 


8-1, 
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6. Peak rock velocity and displacement. Utilizing an 
adaptation of work performed by Cooper (1973), the 
peak horizontal and vertical particle velocities and 
displacements can be represented by equations of the 
form: 


ree Ay We Aly a (8-10a) 
Uy = Ay ~A'y _ 7 (8-10b) 
and 

dy = |By ee vo (8-11a) 
dy = Baa Pi (8-11) 


where Ay, Ay, Ay, Ay, By, By, By, and By are con- 
stants and the other variables have definitions as 
previously defined. 

c. Peak rock stress. The peak horizontal and warical 
stresses are identical to the definitions presented in 
equation 8-5. 

d. Mean peak rock response. The mean peak 
horizontal and vertical accelerations are evaluated by 
using equations 8-1 (variables redefined for rock 
where necessary) with the values presented in 
paragraph 8-3e. The corresponding velocities and 
displacements are obtained from equations 8-10 and 
8-11 where the values of the constants from the data 
presented in table 8-1 should be used. The mean peak 
horizontal and vertical stress may be obtained from 
equation 8-5 where the proper mean values for the. 
properties of particular rocks are considered. 

e. Uncertainty of rock response. Equations 8-6 and 
8-9 may be utilized to compute the uncertainty of 
rock acceleration and stress. The values of 2, and @,, 
presented in paragraph 8-3¢ are still applicable. The 
computation of velocity and displacement uncertain- 
ties produces 


Q2 = Q2 + 4Q2 + E + ail oy 
uy P ; R Ruy 9 
AUC, ZWi/3)]2 “3 , 
+ ee 7 + a (8-12a) 
Vv 
Q2 = Q2 + 402 + 2 + aati" iy 
sg: . Ruy 9 
AiO, ZWU3q2 
+ | -E’ 92 + 92 (8-12b) 
uyR? Z Cy 
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Table 8-1. Mean Parameters for Calculation of Velocity and Displacement Due to Airblast on Rock 


{ 
i a a. © » 06 ZWS oa) ttyl 3 
(o24 2.16 2 3b < 2ywllS 61 eteel/3 
1/3 
a 1193 A,” gsi , o< 203 2 31 tt/ael? 
\u34 10.62 Bie wl 2 er Pepe > 
6 
— (5.03 x 10° sie ds 0 ee zit? 2-siceenel 
5 H j ae in + 
17.82 x 10 -*>8.99 x 10° » dle ae < 61 et/xri/3 
, ‘Ligtiaage | ‘ ee 1/3 1/3 
pg, < p508 x 10° Sal P8099. og, “OZIN 31 -ft/kt 
V 3 3 
"3,35 x 10 13.59 x 10° 4. 38.72 zewi! 3 < 61 £tskrl/3 


_—_— ee 


U.S. Army Corps of Engineers 


yow 
2 —~ 92 4 992 Vv 2 : 
= Dy = a + 925 + 4+ 34 Qy different in character than their vertical counterparts. 
a (1) Typical examples taken from Murrell (1974) 
BiyZWi2 have b ; : 
Vv 92 (8-13a) ave been normalized and are presented at various 
R3dy Z ranges for the vertical velocity in figure 8-1 and for 


Q2 Q2 003 + [4 + SEP op 
Dy ae oe Rid, ) 


2 | me 
Q2 (8-13b) 


BuZWw)2 
+ ER | 
R°dy 
where 2, = 0.36 and Qq = 0.36. 


8-5. Superseismic coil and rock waveforms. 


a, Soil waveforms. According to the data presented by 
Murrell (1974), vertical velocities at a constant range 
maintain a relatively constant waveform but decrease 
in amplitude with depth. Conversely, at a constant 
depth, the vertical velocities undergo distinct changes 
in both waveform and amplitude with variations in 


range. Horizontal velocities qualitatively display the 


same type of behavior although they are markedly 


the horizontal velocity in figure 8-2. The data 
generally apply at all depths in the superseismic 
region. 

(2) Parenthetically, it is noted in figure 8-2 that 
outrunning motion appears at the 340 ft/kt!9 range 
and is characterized by the rolling motion after a time 
of about 1.0 s/kt!/3. The outrunning is well developed 
at the 535 and 680 ft/kt!’3 ranges. The amplitudes 
presented in paragraph 8-3 apply to the motion prior 
to the time of the arrival of the outrunning motion. 
After the outrunning motion occurs, the amplitudes 
of paragraph 8-6 will apply. The reader is encouraged 
to study Murrell (1974) prior to constructing 
waveforms having the generalized waveform charac- 
teristics of figures 8-1 and 8-2 but adjusted to match 
the amplitude statistics presented in paragraphs 8-3 
and 8-6. 
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Figure 8-1. Vertical Velocity Waveforms for Soils at All Depths (Adapted from Murrell, 1974) 
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Figure 8-2. Horizontal Velocity Waveforms for Soils at All Depths (Adapted from Murrell, 1974) 


b. Rock waveforms. Utilizing normalized data from 
Murrell-Carleton (1973), figures 8-3 and 8-4 present 
selections of vertical velocity histories at various 
scaled ranges but at the constant specified scaled 
depths. It is noted that at a given range, the depth 
influences the character of the waveforms (in contrast 
to the data presented in a above in which depth plays 
a minor role in modifying velocity signatures in soil). 
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Conversely, figure 8-5 presents normalized horizon- 
tal velocity histories for a selection of scaled ranges 
and for all depths of interest. The reader should study 
Murrell-Carleton (1973) prior to adapting the data in 
figures 8-3 through 8-5 to suit his own needs. 
Derived signatures should be compatible with the 
peak amplitude data presented in paragraph 8-4 
above for airblast-induced response in rock. 
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Figure 8-3. Vertical Velocity Waveforms for Rock, Scaled Depth of 4.3 fl/kt!/? (1.31 m/ke!’4) (Adapted from Murrell-Carleton, 1973) 
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Figure 8-4. Vertical Velocity Waveforms for Rock, Scaled Depths Below 38.8 fi/kt!”” (11.83 mkt!) (Adapted from Murrell-Carleton, 1973) 
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Figure 8-5. Horizontal Velocity Waveforms for Rock for All Depths (Adapted from Murrett-Carleton, 1973) (1 of 2) 
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Figure 8-5. Horizontal Velocity Waveforms for Rock for All Depths (Adapted from Murrell-Carleton, 1973) (2 of 2) 
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8-6. Outrunning peak ground response. 


a. Peak ground acceleration. Peak vertical and 
horizontal outrunning accelerations (adapted from 
Sauer, 1964) are given by: 


k,2" 


ay = ——(R/W1Sym (8-14a) 
i) 

ay = ay (8-—14b) 

where 


ay = Vertical acceleration, g’s 
ay; = Horizontal acceleration, g’s 
R = Range, ft 
W = Yield, kt 
= Depth, ft 
= Seismic velocity, ft/s 
, 0, m are constants 
b. Peak ground velocity. Equations for the calcula- 
tion of peak vertical and horizontal oscillating outrun- 
ning particle velocities designated as u,, and ap- 
plicable at all near-surface depths of interest in design 
can be represented by the equation 
m = CRW (8-15) 
where 


u.. = Peak oscillatory outrunning velocity, 


ft/s 
c. Peak ground displacements. Peak vertical and 
horizontal out-running displacements for all near- 
surface depths (adapted from Sauer, 1964) are given 
by 


m 


dy Wis _ D(R/W1/3)s 


(8-16a) 
where 
dy = Vertical displacement, in. 
64 = Horizontal displacement, in. 


D, s are constants 
d. Mean peak ground response. Equations 8-14 and 
8-16 may be evaluated for their mean values from 
the quantities: 


3x 1014 -0.44 
3.75 x 1015 -1.35 
1.5 x 1010 -0.44 
1.88 x 1011 -1.35 
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-3.5, 150 < R/W13 < 800 ,1<Z< 16 


tees -3.5, 150 < R/W!3 < 800 , Z> 16 
-2.0, 800 < R/W18 < 3000, 1<Z<16 
-2.0, 800 < R/W!8 < 3000, Z>16 

D = 6400 

s = -1.98 


The mean peak values for the parameters in equa- 
tion 8-15 are presented in table 8-2 (Higgins, 1975). 
For convenience, the data in table 8-2 are cross 
referenced to the appropriate figure (Higgins, 1975). 
The ratio of the initial peak, u,, to the oscillatory 
peak, u,,, is presented in figure 8-6 as a function of 
scaled range. 

e. Uncertainty of oak acceleration. The uncertain- 
ties of equation 8-14 for peak outrunning accelera- 
tions are: 


2 
2 —~ 02 2 202 202 4 M62 (e_ 

Sire ee ee.) eg Veen e aw 17a) 
Quy = Say (8-17b) 


where Q, = 0.8. 

f Uncertainty of peak velocity. The uncertainty 
associated with the amplitude of the oscillatory com- 
ponent of the vertical outrunning ground shock or for 
the maximum horizontal motion for contact or near- 
surface bursts is 


2 
Q2 = Q2 + 7292 + D2 
Um g R 9 W 


(8-18) 


where Q, is the uncertainty of equation 8-15. The 


- values oF 8, are presented in table 8-2 in which the 


data scatter in (Higgins, 1975) is assumed to be log- 
normally distributed. The uncertainty, 2,,,, of the in- 
itial peak vertical velocity u, (depicted in figure 8-7) 
is given by 


2 
2 — o2 202 4 Fi o92 = 
on, Qe + ron + aw (8-19) 


where the values of 2) are presented in table 8-3 
from data presented in (Higgins, 1975). 

g. Uncertainty of peak displacement. The uncertain- 
ties of the displacements in equations 8-16 are 


ple 
Q2 = Q2 + s2Q2 + Us! Q2 (8-20a) 
Vv 9 
Qyer = 24, (8-20b) 
where 9), = 0.6. 
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Figure 8-6. Ratios of Initial to Oscillatory Peak Particle Velocities 
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Figure 8-7. Characteristic Vertical Quirunning Waveforms (Adapted from Higgins, 1975) 
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Table 8-3. Uncertainties of the Initial Peak Velocity 


Geology 


Wet clay over shale 


Dry clay over shale 


Sand over sandstone 
Tonalite 


Clay 


Uncertainty 


“h 
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8-7. Outrunning waveforms. 


a. Velocity. Outrunning is a complicated phenomenon 
that is not easily described by simple waveforms. 
Nevertheless, over the years, certain waveform 
characteristics have been identified that describe the 
principal elements of the velocity signature. The 
reader will be required to adapt the signatures recom- 
mended below to the amplitudes presented in 
paragraph 8-6. . 

b. Vertical. Figure 8-7(a) presents the recom- 
mended waveform for vertical outrunning motion in 
layered soft rock and soil sites. Figure 8-7(b) shows 
the recommended waveform for homogeneous sites. 

c. Horizontal. Horizontal velocity waveforms vary 
continuously with range, and, near the crater, with 
depth as well. Composite waveforms in figure 8-7 are 
more difficult to develop. Figure 8-8 presents typical 
waveforms for wet and dry clays over shale, and sands 
over sandstones. These data are presented for three 
crater-volume scale ranges based on apparent crater 
volumes. Refer to (Higgins, 1975) for definition of time 
scales. The reader may also obtain some guidance 
from paragraph 8-5a on how outrunning develops. 
Consult chapter 12 for the computation of crater 
volumes. 


8-8. Height of burst effects. 


a. Correction factors. Detonation at positive heights 
of burst will alter overpressure and initial coupling 
of energy to the ground, hence the airblast-induced 
ground shock in the superseismic region. To a first 
approximation, HOB-corrected accelerations and 
velocities can be obtained by applying HOB peak 
overpressure correction factors; displacements can be 
modified by applying the HOB-correction factors for 
the overpressure positive-phase impulse. 

b. Velocity, acceleration, and stress. The velocity, 
acceleration, and stress amplitudes characterized in 


equations 8-1, 8-2, 8-5, and 8-10 must be adjusted 
for HOB effects. At.any scaled range in these equa- 
tions, the ratio of the overpressure at the appropriate 
scaled HOB to the overpressure for a contact burst 
from chapter 7 establishes a correction factor used to 
adjust the peak velocity, acceleration, and stress. 
Thus, 


P 
4HOB) = so Ao) eee) 
solo) 
where 
ayop) = Acceleration for a particular 
HOB and range 
P..HOB) = Overpressure at that same HOB 
and range, from chapter 7 
Ps) = Overpressure at that same range 
for a contact burst, from chapter 7 
Ao) = Acceleration at that same range 
for a contact burst, from equation 
8-1. 
and 
P. 
UHOB) = oa Gs (8-22) 
sa(o) 
P 
%HOB) = ee (8-23) 
so(o) 
where 


UWHOB) COB) = Velocity and stress for a par- 
ticular HOB 
Yep» %) = Velocity and stress at that 
same range for a contact 
burst using equations 8-2, 
8-5, and 8-10 
c. Displacement. The peak displacement data of 
equations 8-3 and 8-11 also must be adjusted for 
HOB effects. At any scaled range in these equations, 


8-15 
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the ratio of the impulse at the appropriate scaled HOB 
to the impulse for a contact burst establishes a factor 
that is to be used to adjust the peak displacement. 
Thus, _ 


Ios) 
deop) = = doy (8-24) 
{0} : 
where ; . 
dajop) = Displacement at a particular HOB 
and range 
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Impulse at that same HOB and 

range, from chapter 7 

= Impulse at that same ‘range for a 
contact. burst, from chapter 7 

dj.) = Displacement at that same range 

for a contact burst, from equations 

8-3 and 8-11 


Taio) = 


The ratio of Iqopy1() as a function of range and HOB 
is presented in table 8-4. 
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Figure 8-8. Normalized Horizontal Velocities for Clay over Shale and Sand aver Sandstone Sites (Adapted from Higgins, 1975) 


Table 8-4. Impulse Lost Ho) for Displacement Computation — 


Scaled Range Hop ¢t/kt!/> (mkt ?/4) 


ft/kti/3 (m/e?! 55 200(61) - 400(122) 600(188) 


Near-Ideal Surface 
(91) 
122) 
(183) 
(244) 
(305) 
(381) 


(457) 


(610) 1.1 _ ‘ 
| Nonideal Surface 


(91) 


(122) 


(183) 


(244) 
(305) 
(381) 
(457) 
(610) 
(762) 
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d. Uncertainty. For all practical purposes, the 
uncertainties associated with acceleration, velocity, 
and displacement are equivalent to those established 
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for contact bursts, as defined in paragraphs 8-3f g, 
A for soils, and paragraph 8-4e for rock. 
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CHAPTER 9 


9-1. Gross phenomenology. | 


a. Direct coupling of weapon energy from near- 
surface, contact, and shallow underground bursts 
drives shock waves through the ground and leads to 
crater formation. The relative importance of airblast- 
induced signals, unloading signals from the ground 
surface, and signals resulting from crater formation 
Ff depends on the location of the facility relative to 
ground zero (GZ). As shown in figure 9-1, the region 
between 45 deg (0.79 rad) and 30 deg (0.52 rad) is in- 
fluenced by composite signals of both airblast-induced 
- and crater-induced shock waves. In the central con- 
ical region below depression angles of ~45 deg. crater- 
induced signals are dominant and are the subject of 
this chapter. 
b. The basic variables affecting the intensity and 
duration of the crater-induced ground shock are: 
~ Weapon design 
- Weapon yield 
- Burst elevation or burial 


8@™~ 30 DEG 


CRATER- 
INDUCED 
REGION 


\ 
U.S. Army Corps of Engineers 


CRATER-INDUCED GROUND SHOCK 


- Properties of the site media (rock, soil) 

- Slant range to the facility (or any other point of 

interest) 

c. Crater-induced ground shock can be predicted 
from (1) mathematical simulation, (2) scaling and ex- 
trapolation of near-surface burst data, and (3) scaling 
and extrapolation of fully contained burst data. 
Method 3 currently offers the most credible basis for 
estimating the crater-induced ground shock environ- 
ment and its uncertainty. In this method, under- 
ground nuclear ground-shock environments from 
near-surface (contact or shallow-buried) nuclear 
bursts are estimated from the ground-shock data of 
fully contained bursts. The regions for which these 
data apply include the crater-induced region in figure 
9-1 and, with some reservations, the composite region 
as well. 

d, Of importance to design of facilities are ground- 
shock intensity and time histories and their pulse 
shape. These features are delineated in paragraphs 
9-2 and 9-3. 


~4S DEG 


AI RBLAST-INDUCED REGION 


COMPOSITE 
a 


Figure 9-1, ‘Regions of Ground Shock 
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9-2. Ground-shock intensity. 


a. Decay curves for peak acceleration, velocity, and 
displacement are presented in the general form 
(Perret-Bass, 1975). 


1/3 RM (9-1a) 
aW, = C, wis : 
. 3 ah os 

wis 
R Nn; : . 

13 = 9-1 
dW} C3 Gr 5| (9-1c) 
where ; : 

a,u,dé = Peak acceleration (g), particle 


velocity (ft/s), displacement (in) 


C,n © = Site geology parameters (see table 
9-1) 

R = Slant range from GZ (ft) 

W, . = Effective weapon yield (kt) 


6. In using equations 9-1, which pertain to fully 
contained bursts, a weapon-yield correction is 
necessary to convert to the ground-shock response for 
near-surface bursts. This correction utilizes the 
equation _ : 


W ; : 2 3 
—= : $s 4. B a yas 
W wis w2s Ww 


zt 


+ 
w43 


| a 1/3 
+é | ys Bg < 20 fet 


(9-2) 


where W, and W are the effective and true yields in - 


kt; z is depth (z > 0) or height of burst in ft; a, B, y, 
and é are coefficients of HOB/DOB influence and ee 


accounts for scaling of deep-buried bursts to contact . 
bursts as a function of site material properties and 


the ground motion parameter of interest. 

c. The mean peak of the crater-induced ground 
shock environment is calculated from equations 9-1 
and 9-2, using values of the site-dependent coeffi- 
cients and exponents presented in table 9-1. The 


9-2 


values of the coefficients in equation 9-2 are 
presented in tables 9-2 and 9-3. 

d. The uncertainties associated with the crater- 
induced ground shock for deeply buried systems are 


(+n, 
2=— o2 2 a" 2 2 02 es 
an a oma 5 a QF (9-3a) 
Q2 = Q2 + 92 Mat g2 2 Q2 9-3b 
uu poe Ca 9 oy. t Bu R (9- ) 
(l-nj)? 
2— 02 2 4 3" 92 202 & 
Qe= OF + QE + QW +N Qe (9-3c) 


where {, are the uncertainties of the functional forms 
of the decay laws, 2¢ are the uncertainties of the coef- 
ficients, and so forth. Assume 2, = 0.05 and compute 
Qo, from equation 2-18 where Ly and L, are the up- 
per and lower bounds on the data at any given slant 
range in figures 9-2-through 9-5. Similarly, compute 
Qo and 2, from figures 9-6 through 9-9, and 9-10 
through 9- 43, respectively. Assume n = 1.65. 

e. The uncertainties associated with the yield cor- 
rection presented in equation 9-2 are 


Q2, = O2 + all 1 + 2a2/3W13 + B22/3W23 
e Ww: : 
2 
= eearawse| Q2 
+ 22] a2w4/3 + g272w2/S + 224 


2 
+ eaow-20) Q2. 
+ Z| awe + 282W3 + B8yz2 
(9-4) 


wale 
+ se28/Wi8| 224 + 22 
z "m 


where @, is the uncertainty of equation 9-2, Qy is 
the uncertainty in the yield, 2, are the uncertain- 
ties of ie coefficients, and so forth. Assume 
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Table 9-1. Mean Values of Parameters for Caiculations of Direct-Induced Ground Motion (equation 9-1) 


“Al luvium 66-262 
;Alluvium | 197-1148 

- aWe!/3(g-,1/3 Dry. tuff 131-492 
- Dry tuff (328-1641 
Eq. (9-1a) Wet SEE. - 98-1969 

. Hard rock ° 131-656 
295-7218 


Hard rock 


Mluviun 98-492 

Alluvium 328-1148 
Steere) ‘Dry tuff 131-1641 
am Wet tuff 98-1969 

Hard rock 131-7218 


. 131-492 

3 in/krt/3y Alluvium 328-1148 

Dry, tuff - 328-1641 

Eq. (9-1¢) ia tuff 164-1969 
Hard rock 


U.S. Army Corps of Engineers 


ee 
1/3 


* (100-500) 


230-7218 — 


(20-80) 
(60-350) 
(40-150) 


(30-600) 
(40-200) 
(20-2200) 


(30-150) 
(100-350) 
(40-500) 
(30-600) 
(40-2200) 


(40-150) 
(100-350) 
(100-500) 
(50-600) 
(70-2200) 
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Table 9-2. Factors in Calculating HOB/DOB Effective Yield fequa- 
tion 9-2) 


Coefficient Value 


3.39 x 10 


~3.33 x 10 
2.36 x 10 


29 
-4 
Pi 
-1.19 x for? 


U.S. Army Corps of Engineers 


Table 9-3. Factors [, to Account for Geological Properties (equation 9-2) 


Acceleration 
Geology and Velocity Displacement 


Alluvium 
Tuff 


Hard rock 


U.S. Army Corps of Engineers 
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Figure 9-2. Attenuation of Scaled Acceleration—Dry Alluvium 
(Perret-Bass, 1975) 
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Figure 9-3. Attenuation of Scaled Acceleration—Dry Tuff (Perret: 
Bass, 1975) 
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Figure 9-4. Attenuation of Scaled Acceleration—Wet Tuff (Perret- 
Bass, 1975) 
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Figure 9-5. Attenuation of Scaled Acceleration—Hard Rock (Perret 
Bass, 1975) 
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Figure 9-6. Attenuation of Particle Velocity—Dry Alluvium (Perret-Bass, 1975) 
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Figure 9-7. Attenuation of Particle Velocity—Dry Tuff (Perret-Bass, 1975) 
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Figure 9-8. Attenuation of Particle Velocity—Wet Tuff (Perret-Bass, 1 975) 
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Figure 9-9. Attenuation of Particle Velocity—Hard Rock (Perret-Bass, 1975) 
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Figure 9-10. Attenuation of Scaled Displacement—Dry Alluvium 
(Perret-Bass, 1975) 
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Figure 9-11. Attenuation of Scaled Displacement—Dry Tuff (Perret-Bass, 1975) 
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Figure 9-12, Attenuation of Scaled Displacement—Wet Tuff (Perret-Bass, 1975) 
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Figure 9-13. Attenuation of Scaled Displacement—Hard Rock (Perret-Bass, 1975) 
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9-3. Pulse Shapes 


a. Ground-shock histories are required for the 
analysis of some aspects of system survivability. For 
underground locations that are below a declination 
angle of 30 deg, it is assumed that particle velocity 
can be represented by simple triangles, as defined in 
figure 9-14a. The corresponding acceleration and 
displacement plots are shown in figure 9-145 and 
9-14c. The symbols in the figure are defined as 
follows: 


a = Peak acceleration from equation 9-1a, g 
d = Peak displacement from equation 9-le, ft 
te = 2t,, sec 
t. = 302 a = rise time, sec 
t, 
t_ 


= Positive phase duration = 2d/u, sec 
Negative phase time = 4t,, sec 

u Peak velocity from equation 9-1b, ft/sec 
u, = Peak negative velocity = u/m, ft/sec 


ll 


_R 1/3 
wiss 261 ft/kt 


R 


V3 
wis? 261 ft/kt 


(9-5) 
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6. The mean pulse shape is derived by utilizing the 
above relationships and the data in table 9-1. Also 
utilize 
“a = 0.694 

b = -7.05 x 10-3 

c = 3.15 x 10-5 

c. The uncertainty of a, u, and d are presented in 
equation 9-3 and in conjunction with equation 9-4. 
In addition, assume that 


2 — o2 2 
or = 22 + 92 
Q2 = Q2 
te th 
Q2 = Q2 + Q2 
ty u 
2 ~ 02 
a = 
Q2 = Q2 + Q2 
Uy, u 
al | 
= 2 2 
Qin 0.2 + QF + 5 ew, (9-6) 
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Figure 9-14. Pulse Shape Parameters 
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CHAPTER 10 
CRATER 


10-1. Gross phenomenology. 


a. A crater is formed when a nuclear weapon is 
detonated below ~ 35 ft/ktl/3 (~ 10.7 m/kt!’8),. The 
significant features of the crater include the apparent 
crater (i.e., the visible contour), the true crater as 
defined by the limits of dissociation of the natural 
material, ejecta that falls back into the true crater, 
a rupture zone in which cracking and crushing occurs, 
a plastic zone (or displacement zone in rock) in which 
significant permanent deformation occurs, and the 
ejecta that, along with ground deformation, form the 
crater lip. These zones are depicted in figure 10-1. 

b. The basic variables affecting the cratering 
phenomenon are: 

- Weapon design 

- Weapon yield 

- HOB/DOB 

- Site condition (stratification, water-table eleva- 

tion, etc.) 

- Material properties (principally shear strength) 
The greater the yield, the greater will be the DOB 
(or the lesser the HOB), or the weaker the site 
medium, the larger the crater volume. The crater will 
be flatter the higher the burst or the weaker the site 
medium. Crater size increases with DOB up to the 
optimum cratering DOB (=131 ft/kt™/3, or 40 m/kt9), 
Further increases of DOB reduce crater size. 

c. Currently, crater dimensions are estimated by 
scaling from the cratering data obtained during the 
atmospheric nuclear tests (perhaps a score of crater- 
ing events) and from the large number of atmospheric 


HE tests conducted with charges weighing fromafew . 


grams to 500 tons. The credibility of mathematically 
simulated craters has not been generally established; 
the prediction methods described in this volume do 
not include these data. Moreover, fairly recent data 
from the AEWES, DNA-PR-0012, Phase 2 and DNA- 
PR-0028, Phase 3 (1976, 1977), which included crater- 


ing experiments in soils with a high water content, | 


have not been included in this manual since the 
results were reported after the manual was prepared. 


d. The volume of material excavated by a nuclear | 


blast is of little direct concern to designers except as 
it relates to the formation of ejecta material which 
produces an impact phenomenon and buildup of ejec- 
ta on the ground surface beyond the crater lip. It also 
contributes to the formation of the dust cloud. These 
subjects are discussed in Chapters 11, 12, and 13. A 


detailed derivation of crater volume is presented in 
Chapter 12. A system cannot survive if its original 
position places it within the crater proper or even 


‘ somewhat beyond the crater. Moreover, only very sim- 


ple, compact, inherently rugged facilities can be ex- 
pected to survive in the contiguous rupture zone. The 
crater (crater proper and crater rupture zone) pro- 
duced by a single, relatively small-yield weapon poses 
a threat to a near-surface facility. The crater produc- 
ed by a single super-yield weapon, a laydown pattern 
of weapons, or a sequential detonation of weapons 
delivered to a single point poses a threat to a deep- 
underground facility. The facility designer will thus 
need to know the dimensions and uncertainties of the 
true crater, the rupture zone, and the plastic zone 
under the three patterns of attack. 


10-2. Single burst. 


a. True crater. The true crater is estimated from the 
apparent crater dimensions shown in figures 10-2 
through 10-5 for a variety of soils and rocks. The 
generalized equations for the apparent crater dimen- 
sions are 


R, = a'W!3 + p'z + y'23W-23 : 


Wiis 


< 16.4 


ft, 


2 1/3 3w -2/3 
D,= aW'? + Bz + yz°W 


and 


R= A’W43+ Biz + C'22W-18 


ft, 16.4<—— < 100 
D, = AW!3 + Bz + Cz2W-1/3 we 


;, (10-1) 
and the true crater dimensions are 
Ri = $R, 
D, = $D, (10-2) 


where R is radius, D is depth, z is HOB/DOB (all in 
feet, positive for DOB and negative for HOB), and W 
is weapon yield (in kilotons). The quantities a, f, y, 
¢, A, B, and C are defined below. 

(1) The mean true crater radius and depth are 
calculated from equations 10-1 and 10-2 using the 
values presented in table 10-1. 
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Figure 10-2. Apparent Crater Dimension for Dry Soil or Soft Rock (<10% Water Content), (DNA, 1972) 
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Figure 10-3. Apparent Crater Dimensions for Wet Soil or Soft Rock (>10% Moisture Content), (DNA, 1972) 
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Figure 10-4. Apparent Crater Dimensions for Dry Hard Rock (<3% Moisture Content), (DNA, 1972) 
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Figure 10-5. Apparent Crater Dimensions for Wet Hard Rock.(>3% Moisture Content), (DNA, 1972) 
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(2) The uncertainties associated with the true ' 


crater dimensions are 


Pee ae oe 
OF = E+ BEF QE 


$ 
2 02 2 2 = 
Qh Qe + Qe + 2, (10-3) 
where 
2 —~ 024 LL), 22/3 4 912-2 4 y' 2-6 y7-4/3| 02 
a aes IE a! 2Wele + Br4z4 4 y' 429W |o2. 
+ 22 (p" + By‘z2W-2/8)2 @2 
2/3 
W a’ = 2y'z3W-1)? 2% 5 
Zz 
—, < 16.4 
ws (10-4.a) 
and 
Q2 = 924 LI) a2w24 + B’222 + C’2z4w-2/3|Q2 
Ra £ RZ Co 
a 


+ 22 (B’ + 2C'zW-V/3)2 Q2 


2/3 
+ War crt w-28)? ot 


Zz 
16.4 < —— < 100 
wis 


(10-4b) 


In these equations Q- is the uncertainty associated 
with the functional form of R,, and Lo, is the uncer- 
tainty associated with the constants in the equations. 
Assume that 2; = 0.36, Q2 + Q? = 0.01; assume that 
Qq. = 0.07 for reasonably homogeneous sites and 
0.12 for layered sites or for sites where a water table 
intersects the crater. These uncertainties were ob- 
tained from equation 2-17 where n = 2 and R = 1.35 
for the former case (homogeneous) and R = 1.67 for 
the latter (layered). The calculation of Qe parallels 
equation 10-4 where the notation R, is replaced by 
D.. 
b. Rupture zone. The rupture zone is also estimated 
from the apparent crater, and its dimensions are 


Rp = 6oR, 


Dp = 6'¢D, (10-5) 
Values of 6, 9’, and ¢ are defined below. 

(1) The mean rupture zone radius and depth are 
calculated from the values presented in table 10-1 
and equations 10-1, 10-2, and 10-5. 
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(2) The uncertainties associated with the rupture 
zone are 


2 = 2 2 2 2 
OF OT Oe aOR: OF 
2 = 92 + Q2 4 92 + 9? ts 
oF, = + : a (10-6) 


where Qg, 2g, and 24 are the uncertainties of the 
coefficients in equation 10-5, and Qe, and Qp, are 
defined in a (2) above. Assume that gi + Q2 a Q2 = = 
Q2 + 2 + QP = 0.05. : 
c. Plastic or displacement zane. The plastic zone (for 


soils) and the displacement zone (for rocks) are 
estimated from the apparent crater, and their dimen- 


sions are 
R, = ToR, 
Ds = T'oD, (10-7) 


Values of [ and I' are given below. 

(1) The mean plastic or displacement zone radius 
and depth are calculated from the values presented 
in table 10-1 and equations 10-1, 10-2, and 10-7. 

(2) The uncertainties associated with the plastic 
or displacement zone are 


2 aon 2 2 2 L 2 

ep Qe + Qe + Oe QR, 

Q2 = 92 2 2 2 - 
Dp g + &, + B+ OF (10-8) 


where 9, &-, and 2, are the uncertainties in equa- 
tion 10-7, and QR, and 2p, are defined in a (2) above. 


Assume that Qe ae 2 = = 22 + Q2, +4 Q} = = 0.05. 


10-3. Multiple, pattern Bure: 


Close spacing of charges (less than 1.4 single-charge 
crater radii) enhance both radius and depth of crater- 
ing. Figure 10-6 shows the enhancement of single- 
charge crater dimensions obtained by reducing the 
spacing between charges in a row-charge at optimum 
DOB. The enhancement is proportional to 1//s where 
s is the charge spacing in ft. As the charge spacing 
decreases, the optimum DOB for each charge must 
be increased by the enhancement factor in figure 10-6 
in order to achieve these results. Since the enemy may 
not be able to achieve penetration to the optimum 
DOB, figure 10-6 is an upper bound. Nevertheless, 
it is recommended that the enhanced crater dimen- 
sions be taken as 


Ri = ER, DY = E'D, 
Ry = ER, | Di = BD, 
R, = ER, Di =D, (10-9) 
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where the enhancement factors are of the form 


E = AsR,)¥? (10-10a) 


Similar equations apply to the depth dimensions, i.e., 


E' = A(s/D,)? (10-10b) 
a. Use the mean value of the data range presented 
in figure 10-6 to determine the mean enhancement 
factor. This occurs when A = 1.2 in equation 10-10. 
b. The uncertainties associated with equations 
10-9 and 10-10 are 


QF = 92 + 02 + Lo? + x2 

Qh iS Qe + 92 

Qn = Qe. re 

OF. = oF + O2 + OF (10-11a) 


where 8, is the uncertainty of equations 10-9, and 
@, is the uncertainty of the constant in equation 
10-10. Assume @,= 0.04 and calculate 2, from 
figure 10-6 and equation 2-17 where n = 2 and Ly 
and L, are the upper and lower bounds on the data 
at any selected burst spacing. Similarly, the uncer- 
tainties of the crater depths are calculated from 


QF. = Q + 92 + fQ2 4+ 2, 

2. = Qh. + Q° 

2 = WB + 2 

2 = 28, + 22, ap Q (10-11b) 


10-4. Repeated attack (nail driving). 


a. Repeated attack, sometimes called “nail driving,” 
is the successive firing of shots in the center of craters 


10-10 


formed by preceding shots to enhance cratering. 
Bounds of the true crater data for tests in rock are 
shown in figure 10-7. These data are represented by 
equations of the form 


E” = A+BN (10-12) 
and the enhanced aedheus es 

R;' = E”"R, = E”$R, 

Dy = E"D, = E"$D, (10-13) 


where E” is the increase in true crater, dimensions 
and N is the number of successive shots. The 
parameters A and B are defined ‘in & below. 

b. The mean increase in true crater dimensions can 
be obtained from figure 10-7 or from equations 10-12 
and 10-13, with the constants defined in table 10-2. 

c. The uncertainties of the true crater dimensions 
from the nail-driving attack are 
Q2 = Q24 (AP 4 pana{ Be Q 

Rt f R? 


D,\2 2B2 
Q2,, = 2 + (A2 + B2N2){ <b)" 02 + O22 + 92 
Dt f D t 
(10-14) 


where , is the uncertainty of the form of equation 
10-13, cpand Lon are the uncertainties of the nail- 
driving enhancement factors applied to the radius and 
diameter, respectively. Assume Q¢ = 0.2 and calculate 
Qc, and Qc, from equation 2-17 where n = 2 and 
L, and L, are the upper and lower bounds on the 
enhancement factors in figure 10-7 for any selected 
number of shots. 


10-5. Crater volume. 


Crater volumes in homogeneous and in layered 
geologies are calculated in chapter 12, where the 
depth of burial by crater ejecta is defined. 
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Table 10-2. Parameters for Increase of True Crater Dimensions 
Due to Nail Driving 


Rock 


Radius 


Depth 


Radius 


Depth 


U.S. Army Corps of Engineers 


TM 5-858-2 
CHAPTER 11 
EJECTA/DEBRIS IMPACT 


11-1. Gross phenomenology. 


a. Material is thrown out of the crater region as a 
part of the unloading process following the initial 
strong shock driven into the ground by the explosion. 
Ejected particles produced by a near-surface burst ex- 
tend over a range from the crater lip to remote points. 
Discrete particles are thus part of the crater forma- 
tion, extending through the continuous ejecta 
coverage, on out into the discontinuous region. The 
resulting projectiles of debris create the ejecta/debris 
impact hazard and may pose a substantial threat to 
exposed structures. (Projectiles may also result from 
fragmented structures and trees between the point 
of interest and ground zero, GZ. These are of second- 
ary importance and are not treated here.) Ejecta im- 
pact encompasses the problems of estimating the 
initial mass of ejecta and the size of its consequent 
particles, the number of impacting particles, their in- 
coming velocity and angle, and their distribution. 
Designers must also consider the damage potential 
of ejecta impact in its timing relative to other nuclear 
weapon effects. 

b. The basic variables affecting the number, size, 
velocity, and angle of approach of the impacting 
ejecta/debris projectiles at target points are: 

- Range to GZ 

- Weapon yield 

- Weapon design 

- Height of burst 

- Geology 

- Target elevation 


11-2. Ejected mass and particle size. 


a. The effects of ejecta impact are determined from 
calculations of total ejected mass, number of impact 
particles, mean maximum particle size and number, 
and mean impact velocities and angles. Using the 
data of Post (1974) and Crawford et al. (1974) state 
that the total ejected mass can be calculated from the 
apparent crater mass with the equation 


M,/M, = AW" (11-1) 
where 

M, = Apparent ejecta weight, lb 

M, = Total ejecta weight, lb 

= Yield, kt 
A,n = Constants depending on geologic 
properties 

Since 
M Vv (11-2) 


where 
V, = Apparent crater volume, ft® 
it follows that 


M, = Ay,W®V, (11-3) 


In turn, the maximum particle weight at any range 
(from unpublished work by Post) is determined from 
the total weight to be 


m,, = BM™R¢ , Ib (11-4a) 


where the particle configuration is assumed to be 
spherical, that is, 


m (11-4b) 


m = 7a3,y,/6 


Therefore, the maximum particle diameter is 


an = 1.241 yO eB Awan aRe eva, ft 
(11-5) 


where 


R = Range, ft 
R, = R/V 
B,m,n,a = Constants 


b. The mean value for equations 11-1 through 11-5 
can be obtained from the values of the constants 
presented below: 


tee (hard rock) (ee , We 1kt 
A= 


n= 
0.54 (other geologies) | 0 


,W<lkt 
B = 0.1094 
m = 08 
a = -1.644 
c. The uncertainty associated with equation 11-5 is 
-1)2 ‘ 
Q2 = Q24 (m 1) Q2 (mn) Q2. 
: *P i) 
a m a 
+ —Q2 + |-— - =] 22 (11-6) 
R ‘ :) v 


where 2; = 0.36 and Q, = 0.05. (@{, is defined in 
ch. 12.) 2 
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11-3. Number of impact particles. 


a. It is useful to express the number of particles equal 
to or greater than a specified size that will occur at 
a specified range from GZ. Such an equation is of the 
form 


KsD s-3 _ gs-3 
=ieeee ee (11-7) 
(s-3) an 
where 
I = Number of particles per unit area (ft*) 


greater in size than a, (ft) (an arbitrary 
particle size) but less than a,, (ft) 

s = Constant depending:on geology 

K = by 


ny ; 
= Postattock unit weight of ejecta, lb/ft3 
Equation 11-7 was derived by Post (again, in un- 
published work) based on a review of empirical data; 
it implicitly contains the range parameter from the 
variables D, and a,, where D, is the continuous ejec- 
ta depth, as defined in chapter 12. 
b. The mean number of impacting particles per unit 
area is obtained from equation 11-7. For a, 2 4 in., 


0.00147 for unconsolidated alluvial soil 
S = < 0.5 for hard foal cohesive soils, or layered 
geologies 


c. An estimate of the uncertainty associated with 
equation 11-7 is given by 


-2.5 _ Ga-2.5 
1 [az*” - Gaze \2 


2=— 02 2 aii ee eee AE La 
Ee eee Bede cue as 
ag m 


Q2 (11-8) 
an . 


where Q, = 0.2, Q, m8 obtained from equation 11-6, 
and 2p, is presented in equation 12-2. 


11-4. Particle impact velocities and angles. 


a. Computational and experimental evidence support 
the adoption of a model of ejecta trajectories based on 
the assumption that particles are ejected at a random 
elevation angle 6 and at a velocity dependent on size, 
and that the trajectory is determined by gravity and 
aerodynamic drag forces. Analysis of the trajectories 
obtained from calculations have led to the simplified 
description of impact velocity reported by Crawford 
et al. (1974): 


1/2 
¥: “of, R } exp(-RR,) , R<R,/2 
sin 26 
1/2 
ae pa , R>R,2 
Cy/ (11-9) 


where 
V; = Impact velocity, ft/sec 
Cp = Drag coefficient of particles 
= Range from GZ, ft 

Yp = Preattack density of ejecta, lb/ft? 

90a_y 
R, = —2 sin 26, ft 

D 

ay = Arbitrary particle. diameter, ft 
6 = Ejection elevation angle, deg 


eZ, B = Constants 

6. At an arbitrary range, the incoming particle size 
may be selected in the range a, < a,, where a, is 
obtained from equation 11-5. The ejection elevation 
angle is assumed to be distributed between 15 deg and 
75 deg, but strongly concentrated near 45 deg, as 
shown in figure 11-1. The drag coefficient is assumed 
to lie between 0.3 and 1.2, with a mean value of 0.6. 
The mean impact velocity is computed from equation 
11-9 where .W = 5.673 and B= 23.08. 

c. The uncertainty in the impact velocity is given 


by 
2 
Q2+ te = BS Q2 4 20 _ 1 eot229 92 
0 


0 


; 
ad Ron 1 (Bl oe 
*y, +] a + (| mep,. © peste 
1 
24 —O2 2 
+ 72 +79 » R>R/2 
(11-10) 


where &- = 0.36 and Qe = M% = 0.3. 


11-5. Impact damage potential. 


a. Impact damage potential can be measured in 
terms of particle kinetic energy (momentum), ex- 
pected target penetration, or expected crater volume. 
For hardened facility applications, use the damage 
measures—penetration or crater volume. Penetration. 
has been examined as a continuing part of armor/ 
armor-piercing projectile development and, more 
recently, as part of the safety program for nuclear 
power reactor development. Available data for pro- 
jectiles striking concrete targets are utilized to ex- 
trapolate the effect of rock fragments striking similar 
targets. 

b. The penetration depth, d, is adapted from ASCE 
(1976): 


0.04 Fp(KNy,)? altvio9 | da, <2 


2 x 10-6 FpKNy,ai2V;18 , d/a, > 2 (11-11) 


FRACTION OF TOTAL MASS 
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@, EJECTION ANGLE DEG 
Figure 11-1. Ejection Angle Distribution for High-Yield Surface Bursts (Seebaugh, 1975) 


where 
Fp = Particle deformability factor 
_ §0.7 for hard rock 
Hee for soft rock or alluvium 


Pate 180 
K Concrete penetrability index VE 
f, = Ultimate compressive strength of con- 
* crete, psi 
N = Particle nose shape factor = 1.0 
a, ,= Arbitrary particle diameter < a,,, in. 
V,; = from equation 11-9 based on particle 


size a,, ft/sec 
Yen Preattack particle density, lb/in.? 
The estimate of penetration crater volume is based 
on an assumed conical crater depth d and radius 2.5d 


@ V crater = 6.5 d3 (11-12) 


c. The uncertainty associated with equation 11-11 
is 


28 1 1 
2 see ,loe loge 2 
yt TMK + Tt TA + 0.894, , 
d/a, <2 
Q2 = Q2+ ; 
df d-ao\?} o5 24092 102 
Fi Pap, +9 + ok a8 
+ 3.24 <4, da, >2 (11-13) 
1 


where 9, = 0.5, Qn, = 0.25, Qy = 32, Qh = 0.15, 
Q., = 0.05, and Qy, is obtained from equation 11-10. The 
quantity Qe is the uncertainty of the unconfined com- 
pressive strength of concrete. 
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CHAPTER 12 
EJECTA/DEBRIS DEPTH 


12-1. Gross phenomenology. 


a. To prevent burial of a facility and interruption of 
its function or to plan active removal of debris, the 
designer must know the expected depth of debris, in 
addition to the expected particle size distribution 
delineated in chapter 11. Deposition of crater ejecta 
as surface debris becomes significant when the height 
‘of burst is less than 35 ft/kt!-3 (10.7 m/ktl9), The 
depth of debris in the region of continuous coverage 
of the ground depends on range from ground zero (GZ) 
and varies (usually unpredictably) with azimuth, 
probably the result of fracture patterns and in- 
homogeneities in the medium. Outside the continuous 
region (greater than five crater radii), the debris con- 
sists of separated discrete particles unlikely to bury 
the facility. The amount of ejecta threatening burial 
of a facility depends on the crater volume and is ad- 
dressed in this chapter for both homogeneous and 
layered geologies. 

b. The basic variables effecting the ejecta thickness 
at a point on the ground are: 

- Range to GZ 

- Weapon yield 

- Weapon design 

- Height of burst or depth of burst 

- Geology 
Ejecta depth is enhanced by large yields, reduced 
ranges, weapon penetration, and certain geological 
features described below. 


12-2. Debris depth. 


a. The depth of crater ejecta debris from a near- 
surface burst is related to the apparent crater volume. 
(See chapter 10 for calculation of apparent depth and 
radius.) In turn, the crater volume is implicitly depen- 
dent on weapon yeild, height of burst, weapon type, 
and geology (Post, 1974). The scaled debris depth as 
a function of scaled range can be represented by an 
equation of the form 


De R R 
yiss =A exp Coie te (p vis ) (12-1) 
a a a 
where 
D, = Average ejecta thickness, ft 
V, = Apparent crater volume, ft? 


R = Range from ground zero, ft 
A, B, a, and # are constants 


b. The mean average ejecta debris depth is 
calculated from equation 12-1 where the constants 
have the values: 


A B a B 
0.75 0.017 -2.3 -0.423 
Equation 12-1 is plotted in figure 12-1 as the median 
value with uncertainty bounds. These data are ap- 
plicable for depths of burst less than 17 ft/kt!/3 
(5 m/kt!/8) and for most of the layered geologies en- 
countered within the continental United States. 
Equation 12-1 should be multiplied by the quantities 
in table 12-1 for the other geologies shown therein. 
The mean apparent crater volumes are calculated 
from the procedures presented in paragraph 12-3 
below. 

c. The uncertainty associated with equation 12-1 is 


Q2 =~ 924 _1 
9D2 


e 


2 
_ppneon sf 98, 


[p.- AaRe*®Va¥9 


2 2 
+ a eae.) + party, Q2 
D2 (12-2) 


e 


where {2-is the uncertainty of the functional form of 
equation 12-1, Qy, is the uncertainty of the apparent 
crater volume, etc. Calculate 2, from equation 2-20 
where L, and Ly are the lower and upper bounds 
from figure 12-1, multiplied by V1/3. The uncertainty 


of V, is obtained from paragraphs a (2) and a (3) 


below. 


12-3. Crater volumes. 


a. Homogeneous geologies. Normalized apparent 
crater volumes are presented for both “high-yield” 
and “low-yield” nuclear devices. These terms from 
Cooper (1971) refer to “hot” devices having a high 
yield-to-mass ratio and to “cold” devices having a low 
yield-to-mass ratio. The relationship of volume to dif- 
ferent ranges of HOB/DOB are also considered for 
both high-yield and low-yield devices. Typical crater- 
ing efficiencies are shown in table 12-2 for half-buried 
HE charges. 
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Figure 12-1. Ejecta Depth Prediction Curves for Nuclear Events, 1 <R/V!" < 12 (Post, 1974) (1 of 2) 
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Figure 12-1. Ejecta Prediction Curves for Nuclear Events, 0 < R/V "9 < 40 (Post, 1974) (2 of 2) 
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Table 12-1. Geology Correction Factors Pz = 


Hard Rock 


Porous Rock 


U.S. Army Corps of Engineers 


Table 12-2. Normal Cratering Efficiencies E., for Half-Buried HE 
Charges (Post, 1974) : ; 


Idealized Geology 


. Wet Soil 


Wet Soft Rock 


Dry Soil 


Dry Soft Rock 


Hard Rock : ‘ 


12-4 


(1) The mean apparent nuclear crater volumes for 
100 < W < 10 kt can be described by the following: 


“High-yield” weapons: 


V,/V, = 0.04 e~ 55.1, H,2 0.002 

V,/V, = 0.05 e~ 166 Hy 0 < H,< 0.002 

V,/V, = 4.07 (0.0033 + D)°77 0 < D.< 0.02 
(12-3) 


“‘Low-yield” weapons: 


V/V, = 4.39 (D, + 0.0216)0-88 
V,/V, = 4.07 (0.0033 + D077 


0< D,<0.055 


0.055 < D,< 0.2 
= (12-4) 


where 
H, = HOB/V,18 


D, = DOB/V,}4 
V, = WE, . . — 
E, = Cratering efficiency in ft®/kt 

W = Yield in kt 
in which HOB and DOB are both expressed in ft. For 
yields greater than 10 kt, the DOB cases of equations 
12-3 and 12-4 should be replaced by the following: 
“High-yield” weapons: 

ou NOUR 

V,/V,, = 4.07 we[0.0083 + D,wes| 0<D,<0.2 
| (12-5) 


“Low-yield’” weapons: 


0.88 
V,/V, = 4-39 wD, Wa + 0.0216] 0<D,<0.055 


0.77 
V/V, = 4.07 we D,wee + 0.0033] 0.055<D,<0. 2 


(12-6) 


where a = —0.3 D,W?.02 
(2) The uncertainties associated with equations 
12-3 and 12-4 are specified as 


(12-7) 


Q2 = Q2. 1 — THs? ge + H%?o2 
3 Vo s’ “HOB 


for the HOB conditions and 


27,2 
yD, [oe Y D: y) 


Q2. = Q2 +(1 -——- ——_=+—_ ¢ 
Va f 3(A + D,) Vo (A D,)? DOB 


(12-8) 
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for the DOB conditions. The values of the constants 
A and y are defined as: : 


Assumed 
Condition Values of Constants 
A y 
“High-yield”” weapons 
HH, 2 0.002 - -55.1 
O0< Hy, < 0.002 - -166.0 
O< D, < 0.02 0.0033 0.77 
“Low-yield” weapons 
O0< D, < 0.055 0.0216 0.88 
0.055 < D, < 0.2 — 


0.0033 0.77 


The quantity 2, is the uncertainty of equations 12-3 
and 12-4, &hyop is the uncertainty of the HOB, etc. 
Assume 9, = 0.2 and let 


QF = oF, + oF (12-9) 


where Qj, = 0.3 based on the use of equation 2-20 
for an uncertainty band of + 50% on the data 
presented in table 12-2. 

' (8) An approximation of the uncertainties 
associated with equations-12-5 and 12-6 is 


yaD,We3 
2 = 02 — 2 
Sy ep Fr + Dw) se Ris 
yD, We i 2 
(A + DW) 
af3 
ares Sa) 2 (12-10) 
3(A + DW) Vo 


where @r is the uncertainty of equations 12-5 and 
12-6 and y is the uncertainty of V,. The constants 
for the various conditions are: 


Assumed Condition Value of Constants 


A y 
“High-yield” bursts: 0.0033 0.77 
“Low-yield” bursts: 
0< D, < 0.055 0.0216 0.88 
0.055 < D, < 0.2 0.0033 0.77 


Assume 2; = 0.2 and use equation 12-9 to compute 
Qy . 
oO 
12=5 
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b. Layered geologies. Apparent crater volumes can 
be determined for layered geologies from the equation 


Va = Vi V3 
————— = l1- —§.4d/V = 
VeeVi exp ( a (12-11) 
where 

d = Layer thickness in ft 


Vy= Crater volume (at the HOB or DOB 
of interest) when d = © 


V_= Crater volume (at the HOB or DOB 
of interest) when d = 0 


Equation 12-11 can be written for a layered site of 
n different materials via the recursive equations: 


(aa ¥alfts-¥4 


=1-exp [-5.44,s[Vn-aa}™| (12-12a) 
(eee ~V-1f{Va-2-Va-a} 


= 1~ exp [5.44 -2fVn-2n-aa}9| (12-12b) 


1nf{vi- Ve, soe 


ie) (12-12n) 


(V2, pea V2 ae 


= 1- exp -5.44 (V2, . es 
Equations 12-12 must be solved sequentially, begin- 
ning with the first equation and working upward 
through n materials. 

(1) For convenience, the mean for equations 
12-12 coupled with the equations and figures dis- 
cussed in previous subparagraphs are presented in 
figures 12-2 through 12-6 for a few interesting 
geologies. In figures 12-2 through 12-4, cratering ef- 
ficiencies can be determined for three different 
geological conditions. Based on these cratering effi- 
ciencies, the mean ejecta-debris depth can be 
calculated from figure 12-5 for airblast pressures of 
600 and 1000 psi and in figure 12-6 for a larger selec- 
tion of overpressures. (Refer to figure 7-10 to deter- 
mine overpressure/range relationship.) 

(2) The uncertainty of ejecta depth for layered 
sites should be calculated from equation 12-2. Equa- 
tion 12-2 in turn requires the uncertainty Qy,- 
Depending on the HOB or DOB and the yield and type 
of weapon involved, Qy, should be determined from 
one of equations 12-7 through 12-10 (omitting equa- 
tion 12-9, which is used to compute the uncertainty 
of Vy). Equations 12-7, 12-8, and 12-10 require a 
value of V,,. It is recommended that V, be calculated 
from the applicable condition in equations 12-3 
through 12-6 by iterating the value of V,, until the 
ratio V,/V, on the left-hand side of the equation is 
compatible with the right-hand side that involves 
scaled HOB’s and DOB’s. 
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Figure 12-2. High-Yield Nuclear Surface-Burst Cratering Efficiencies in Layered Geologies, Soil over Soft Rock (Shale or Sandstone) 
with a Water Table (Post, 1974) 
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Figure 12-3. High-Yield Nuclear Surface-Burst Cratering Efficiencies in Layered Geologies, Soil over Hard Rock (Limestone) with a 
Water Table (Post, 1974) 
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Figure 12-4. High-Yewld Nuclear Surface-Burst Cratering Efficiencies in Layered Geologies, Alluvial Soil over Hard Rock Baser 
with a Water Table ‘Post, 1974; 
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CHAPTER 13 
DUST 


13-1. Gross phenomenology. 


a. Dust is lofted into the air by the deposition of the 
initial radiation on the ground surface, by the crater- 
ing process, by the outward sweep of the airblast 
across the ground surface, and by the inward sweep 
of the afterwinds produced by the rising fireball. The 
lofted dust is quickly partitioned into a low-lying or 
base surge cloud and a high-altitude or ejecta dust 
cloud. 

b. Only recently have dust clouds been accorded at- 
tention by the nuclear weapon effects community. The 
material presented in this chapter is a synthesis of 
the literature on the subject (see, for example, Young, 
1965; Tewes, 1970; Crawford, 1970; Allison, 1971; 
Swatosh, 1971; Powers, 1975). 

c. The basic variables affecting the potential 
destructiveness of the free-field dust environment 
received at a ground surface location of interest here 
are: 

- Weapon yield 

- Burst elevation 

- Dustiness of the ground surface 

- Dustiness of the crater ejecta 

- Prevailing winds 

- Height of temperature inversion (mixing depth) 
Dust effects increase with larger weapon yields, lower 
heights of burst, dustier ground surfaces, calmer 
winds, or lower inversion bases. 

d. Personnel will rapidly become incapacitated if 
they breathe the heavily dust-laden air. Appropriate. 
design measures must be taken either to remove the. 
dust from the air or to allow for functioning in a 
sealed facility. Knowledge of two aspects of the dust 
cloud are required to implement dust removal 
measures or button-up procedures: 

‘ Particle-size distribution vs. time 
Load vs, time 
For calculations the ground-surface dust and the 
ejecta dust are treated separately as though they 
formed two noninteracting clouds. 


13-2. Base-surge dust cloud. 


a. Cloud model. The model of the base-surge dust 
cloud is developed on the assumption that it is instan- 
taneously formed and initially centered at the weapon 
aim-point. The model is configured as a right circular 
cylinder of constant radius R and constant height H 
and in which the initial dust concentration 
throughout the cloud, w,, undergoes no dilution or 
diffusion of the cloud with the surrounding air but 
decreases in concentration due to vertical settlement 
of the dust particles. The cloud will move horizontally 
at the prevailing wind velocity. 


(1) The cloud radius is given by 


R = R,we (13-1) 


ee) 
| 


= Cloud radius, m 
R, = Cloud radius for a 1 Mt burst, m 
W = Weapon yield, Mt 
n = Constant 
(2) The mean radius is obtained from equation 
13-1 where R, = 3000 m and n = 0.4. 
(3) The uncertainty associated with the radius R 
is 
2 = o2 2 202 ie 
QF QF + i +n Oy (13-2) 
where 9, is the uncertainty associated with the func- 
tional form of R, Qp, is the uncertainty parameter 
pasa with R,, etc. Assume 2 = 0.2 and 
= 0.2. 
Rhy The cloud height is given by 
H = H,W™ (13-3) 
where 
H Cloud height, m 
H, = Cloud height for a 1 Mt burst, m 
_ In = Constant 
(5) The mean height is obtained fram equation 
13-3 where H, = 100 meters and m = 0.2. 
(6) The uncertainty associated with the height H 


is 


or = O2 + oF + m?Qe, (13-4) 


Assume Qe = 0.2, Qu, = = 0.2. 

b. Initial dust concentration. Assume that the top 
£ meter of dust cover is uniformly dispersed into the 
cloud. Accordingly, the initia] dust concentration is 


Wo = * (13-5) 
where 
— = Thickness of dust layer on the earth’s 
surface that is dispersed into the cloud, 
m 
@ = Density of the in-situ earth layer, g/m® 
w, = Initial dust concentration, g/m? 
H = Dust cloud height (eq. 13-3) 


(1) The mean initial dust concentration is - 
obtained from equation 13-5 where ~ = 0.5 m or 
whatever depth applies to a particular site, and 
@= 1.0 x 108 g/m’. 

13-1 
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(2) The uncertainty associated with the initial 
dust concentration is ; 
2 = 92 2 2 2 a 
2, Ort Qe + 2° + QE (13-6) 

where Q2 = 0.4, 2, = 0.4, and Q, = 0.2. 
c. Vertical settlement speed. The vertical settlement 
speed of a dust particle is 


oo ; a > 1000 pm 
v{*] . 100<a< 1000 um 
1000 
eS a \3-0 
2] , 10<a¥< 100 um 
100 
0 ; a<10um 
(13-7) 
where 
a = Particle diameter, um 
V,, Vi, = Settling speeds for a 1000 pm and 


100 pm particle, respectively, m/s 
The meaning of V = ~ for a > 1000 pm is that parti- 
cle sizes larger than 1000 um are not lofted. 

(1) The mean vertical settlement speed of a par- 
ticle of diameter a is obtained from equation 13-7 
where V, = 55 m/s and V; = 5.5 m/s. Mean settle- 
ment speed is-plotted in figure 13-1. 


(2) The uncertainty associated with the settle- 
ment speed is ; 


old = 0.13 (13-8) 


d. Wind transport. Assume that the cloud is 
transported as a rigid body by the prevailing wind. 
Assume that the cloud is initially centered above 
ground zero and thereafter moves with the prevail- 
ing wind. Use local wind statistics (contact the 
National Weather Records Center, operated by the 
Environmental Services Administration at Ashville, 
NC, for meteorological data) to determine transport 
speed and direction of the cloud. Track all clouds that 
originate in the general vicinity of the point of in- 


terest to develop probabilities of their passing over _ 


the point of interest. Appropriately account for the 
“contribution”’ of each cloud. 

e. Particle size vs. time. With the passage of time, 
progressively smaller and smaller particles settle out 
of the cloud. In effect, the cloud progressively exhibits 
a finer distribution of dust. At time 


13-2 


0 ; a > 1000 pm 
H 0.2 : . 
ee bib , 100<a< 1000 um 
Z | 
i] 
t, = 
0.2 ' 
mie ,  10<a<100 um 
(fa) 
*\ 100 
oo : a<10um 


(13-9) 


only particles smaller than diameter a are in the 
cloud. 

(1) The mean value t, is obtained from equation 
13-9 where the mean value of H and V are obtained 
from a (5) and c (1) above, respectively. The mean 
value of equation 13-9 is plotted in figure 13-2. 

(2) The uncertainty associated with t, is 


2 — 02 2 2 = 
fa QF + QF + Qs, (18-10) 


where $2}, and @, are calculated from equations 13-4 
and 13-8. Assume 9, = 0.36. 

f Load vs. time. The dust load that is entrained in- 
to the facility is of critical importance to the designers. 


The load is determined by integrating the product of 


the dust concentration and the air-intake-volume rate 
over the elapsed time of interest: 


fe 
L= pla,t)o(t)dtda (13-11) 
ay ie] 
where 
L = Dust load, g 
t = Time, s 
o = Air-intake-volume rate, m®/s 
wy = Dust concentration at an arbitrary 
time t, g/m? 
@,, a9 = Size of smallest and largest particles 
in cloud, m 


For the assumption of a constant air-intake-volume 
rate and of a particle diameter distribution of a, the 


TM 5-858-2 


io) 


V, MEAN SETTLEMENT SPEED, m/s 


0 7 10! 107 10 
a, PARTICLE DIAMETER, um 


U.S. Army Corps of Engineers 


Figure 13-1. Mean Vertical Settlement Speed -.of Dust Particles. 
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Figure 13-2. Mean Settle-Out Time of Dust Particles—Base Surge 
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load attributable to particles of diameters a, < a 
x ao is 


ao ty 
vs f f wla,t)d(t)dtda 
ax to 
+ f. J wla,t)p(tdtda 
a, Jt 
Axx ts 
+ f J w(a,t)M(t)dtda 


a, rt 
‘ a8 a w(a,t)(t)dtda (13-12) 
a, 4 ty 
where 
7 aH 
ax = —, 
Vit 
1/3 
aux = 10 a2) ; 
Vit 
a, = Smallest particle that settles, m 
ty = 7 : 8 
oO 
ty = = a) 
oO 
_ 1000H 
tz = y ,s 
oO 
y = —%  gim4 


» gh 
aH? n(a,/ay) 


(1) The mean load is obtained from equation 
13-12 by substituting the appropriate mean values 
for variables where a, = 1 wm ag = 1000 um. 

(2) The uncertainty associated with the load L is 


2= 92 2 2 2 2 
Op = Bet Qi + Qe + Qe + Me 


where 2 = 0.4, Q, = 0.2, oe = 0.4, and SH is 
calculated from equation 13-4. 


(13-13) 


13-3. Ejecta dust cloud. 


a. Cloud model. The ejecta dust cloud model assumes 
that the cloud instantaneously reaches a height equal 
to the base of the upper atmosphere temperature in- 
version (H’) and remains at that height, and that it 
instantaneously reaches a radius equal to twice the 


me n(2R,)°H' 
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apparent crater radius and thereafter expands 
radially according to the equation 


(=) 
v= vj 
R’ 


(13-14) 


where R, = apparent crater radius, m; R' = cloud 
radius at time t, m; v = cloud radial surge velocity 
at R', m/s; v, = cloud initial radial surge velocity, 
m/s. The cloud is also assumed to be initially centered 
at the weapon aim-point and has a uniform initial 
dust concentration throughout the cloud, y,, which 
undergoes no dilution or diffusion of the cloud with 
the surrounding air, except for the horizontal expan- 
sion stated in equation 13-14. The cloud maintains 
a uniform, horizontal dust concentration at all times 
but the dust concentration decreases due to vertical 
settlement. of the dust particles. The cloud moves 
horizontally at the prevailing wind velocity. 
b. Radius. The cloud radius is 


R’ = (8R3 + 12v,R*t)/3 (13-15) 
a a 
where 
R, = Apparent crater radius, m 
R’ = Cloud radius at time t, m 
t = Time,s 
v, = Cloud initial radial surge velocity, m/s 


(1) The mean cloud radius at time t is obtained 
from equation 13-15 where the mean value of R, is 
defined in chapter 10 and V, = 200 m/s. The cloud 
will stabilize when the mean surge velocity is very 
smal] compared to the mean prevailing wind velocity. 

(2) The uncertainty associated with phe cloud 
radius at time t is 


Q2,= 07+ 92 +92 (13-16) 


a a 

where , is the carelanarh associated with the func- 
tional form of R', Q, is the uncertainty associated 
with v,, etc. Assume o,= = 0.36 and Q, = 0.36. Refer 
to chapter 10 for calculation of Qn 

c. Height. The cloud height is equal to the height 
of the base of the temperature inversion. The mean 
height (based on a one-year interval) and the uncer- 


- tainty parameter of the inversion base can he deter- 


mined from weather statistics. 
d. Initial dust concentration. Assume that a frac- 


| tion of the earth material that comprised the apparent 
’ crater volume is uniformly dispersed into the initial 


ejecta dust cloud. Accordingly, the initial dust con- 
centation is 


, — farfioooVar : (13-17) 
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where 

fy, = Fraction of the ejecta dust which is 

partitioned into the dust cloud 
= 0.5 - HOB/500 , -250 < HOB < 250 

m/Mt}/3 

fioo9 = Fraction of the crater ejecta compris- 
ing 1000 ym diameter particles and 
smaller 

HOB = Scaled height of burst, m/Mt!? 

H’ = Height of the inversion layer above 
the ground surface, m 

V, = Apparent crater volume, m3 


y = In situ earth density, g/m® 

(1) The mean initial dust-concentration is 
obtained from equation 13-17 where the mean values 
of V, and R, are defined in chapter 10. Mean values 
of fjo99 and H' are site specific. 

(2) The uncertainty associated with the initial 
dust-concentration is 


Q2,= Q2 +92 +92 492 +924 92 +92 
#o fp" fui fi000 Va y Ra H 


(13-18) 
where 
HOB/500_—s\? Sh 
Q2 + Q2 jee 2 (18-19) 
fHL fe \0.5 - HOB/500/ HOB 


and in which 2. = 0.36, 2, = 0.17, and Qs, = 0.36. 

e. Vertical settlement speed. Paragraph 13-2c (for 
base-surge dust) is applicable. 

f Wind transport. Paragraph 13-2d is applicable. 
However, remember that the ejecta dust cloud radius 
is not constant but is increasing with time as defined 
in the model a‘above. 

g. Particle size vs. time. Paragraph 13-2e is ap- 
plicable when H is replaced by H. 


13~6 


h. Load vs. time. The form of equations 13-11 and 
13-12 remain unchanged; however, utilize the follow- ; 


ing values of the parameters: 


ae: Se 
Vit 
eet aad ee (2s of 
10 \V,'t 
ty = H’V,, s 
tg = 10HV,s i 
tz = 1000H/Vi,5 | 
fut fiooo Vay : 
a ~ aH’ R? (e + =)" 1n ad/a, 
R, = 12V,13 


In (ap/ay) 


f 0 = 
wee ae (ap/ay) 


(1) The mean value of the ejecta dust cloud is ob- 
tained from equation 13-12 when a, = largest lofted 
ejecta particle = 18 m, a; = 1 ym, ag = 1000 um. 

(2) The uncertainty associated with the ejecta 
dust cloud is 


1 
2°— o2 2 2 2 2 i 
2 Qe ae ee Qe + QQ + 92 (13-20) 


where %¢ = Qh, = 0.36, 2p = 0, 2, = 0.05, and 2y_ 
is obtained from equations 12-9 or 12-10. 
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CHAPTER 14 
FIRESTORM 


14-1, Gross phenomenology. 


a. A firestorm, which can develop from the uniting 
of the many small fires started by the thermal radia- 
tion pulse, represents the most violent type of mass 
fire. The thermal convection column and strong 
ground-level winds blowing inward restrict growth of 
the fire, but intensify the burning. Firewhirls or fire 
cyclones may develop during the intensified burning 
phase. Full development of a firestorm can occur in 
an hour; burnout may not occur for several hours. 

6. A firestorm results in oxygen depletion, eleva- 
tion of ground and air temperatures (several hundred 
K), reduced or enhanced atmospheric pressures, and 
massive air pollution. Appropriate design measures 
must be taken, either to alter the site fuel-density 

characteristics sufficiently to reduce the probability 
of a firestorm to a safe level, or to provide the facility 
with heat shields and semi-closed-cycle ventilation, 
life support, and combustion air systems for the dura- 
tion of the firestorm. The designer must know the den- 
sity of the site fuel, the types of fuel, and the condi- 
tions that will ignite the fire in the thermal radia- 
tion environment (chap. 6). 

c. After a nuclear burst, the likelihood that a 
firestorm will develop is a function of: 

-. Weapon yield 

- Burst elevation 

- Slant range 

- Weather conditions 

- Fuel density and kindling type 
The first three factors primarily affect ignition 
through their influence on incident radiant energy 
(see paragraph 6-2, Radiant Exposure). The other fac- 
tors account for the availability and condition of fuel. 


14-2. Weather conditions. 


Firestorms can occur in the dry season and if wind 
velocities are less than 8 ft/s (2.4 m/s), a calm wind. 
Climatological data of a facility site are required to 
estimate need for protection and level of protection 
against firestorms. 


14-3. Fuel density. 


a. Conditions. Prerequisites for a firestorm are: 

- Gross area aflame > 300 acres 

- Net fueled area + gross area > 25 percent 

_ Firebreaks < 50 ft wide 

~ Fuel density (based on gross area) > 175 tons/acre 

- Sufficient kindling/tinder available so that 10 
cal/em? of incident nuclear radiation will cause 
ignition 


A near-surface burst of > 10 kt will produce 10 
cal/em? (4.184 x 10° J/m2) radiation intensity at suf- 
ficient distances to ignite 300 acres for virtually all 
dry-season preattack conditions. 


b. Spontaneous ignition. Spontaneous ignition of 
kindling occurs when the incident energy density ex- 
ceeds specific values. Thus when 
QA: > Q% (14-1) 
ignitions will occur; if the conditions specified in a 
above are present, the fire will become a firestorm. 
The value of 9, in equation 14-1 is obtained from 
equation 6-16 in chapter 6 of this manual for a time 
of trax a8 defined below. The value of {2p is also deter- 
mined below. 


(1) Ignition thresholds for dry kindling materials 
are expressed with equations of the form (Kerr, 1971) 


K.H,t 
Qo = -K, ech [+-K,}» + Semen (14-2) 
where 
Qo = Ignition threshold radiant exposure, 
cal/em? 
ec = Volumetric heat capacity, cal/°C/ 
em? 
= Thickness of the material, cm 
a = Absorptivity of material 
x = ¥ 4tmax 
L 
a = Thermal diffusibility of the material, 
cm2/s 
H, = Critical irradiance, cal/em2/s 
tmax = Time to maximum temperature 
= 0.032 W295, s/ktl/? 
WwW = Weapon yield, kt 


K,, Ky, Kg, and n are constants 


For moist fuels, equation 14—2 is modified by the fac- 
tor 1+ hVoL so that 


J 


Qo, = (1+ h¥eL)Qy (14-3) 


where h = relative humidity in percent. 
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(2) Mean values for parameters related to wild- 
land fuels are presented in table 14-1 for use in equa- 
tions 14-2 and 14-3. In addition, the following restric- 
tions prevail: , 


-750, 1, 0, 
K, = ¢ -2250, K, =<0, Ky =<0, 
-1522, 0, 2.6, 
aye + <0.25 
n= 1, 025< 1<0.52 
0.4, t 20.52 


The value of H, is generally not a fixed value for all 
fuels but is dependent on the opacity and geometry. 
However, a value of 0.4 cal/em2/s (1.7 x 104 J/m?/s) 
can be used in lieu of specific values measured in the 
laboratory. 


(3) The. uncertainty associated with equation 
14-2 can be reasonably expressed by 


Weel 2 
ot sa oop PION. Joey 2] ny] 2 
Qo of Q? a? rL 


+4+I1L- KaH tax 2 
K,ecI™ 


a KH. |” 
t2 n a _ 3° *c 
onl t max =a 


14-2 


4 mM tmax + K3L%n? 
472 r2 
2K2H*t? 
+ c.max Q : 
(gc)*K2r2n { * (14-4) 
in which [ = ———— We - Kg 


and where Q, is the uncertainty of the functional 
form of equation 14~2 or 14-3, Q,, is the uncertain- 
ty of the constant K,. For moist kindling, and to a 
first approximation, the uncertainty associated with 
equation 14-3 is 


Q20 = Q2 4 


14-5 
2, = 20 (14-5) 


where 2 is taken from equation 14-4 and Q, is the 


uncertainty of the value of the humidity. 

(4) Using existing experimental data, it is not 
practical to evaluate the various uncertainties 
associated with the various parameters in equations 
14-2 and 14-3. Instead, it is recommended that the 
calculated values for Qy and Qp_ be increased by the 
factor 1.15 and that @, = 0.2 be used based on the 
data presented in table 14-2, in which the mean and 
standard deviation (exclusive of the wheat-straw data) 
were used to estimate the bias and scatter between 
calculated and experimental data. Assume 2g, = 0.2. 
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Table 14-1. Physical and Thermal Properties of Wild Land Fuels (Kerr, 1971) 


Thermal 
Diffusivity, 


Fuel Materials 


Grasses and Other Fine Fuels 


Horsehair lichen (Aleetoria jubata) 


Cheatgrass seedpod (Bromue 
‘tectomen) 


Wiregrass (Aristrida stricta) 
Desert stipa (Stzpa spp.) 
Cheatgrass leaves (Bromus tectorwn) 


Harding grass (Phalaris tuberosa 
var stenoptera) 


Sedge (Carex geyert) 
Wheat straw (fritiewn spp.) 
Cheatgrass stalks (Bromus tectorun) 


Hardwood and Shrub Leaves 


Beech (Fagus spp.) 
Chestnut oak (Q. montana) 


Scarlet oak (@. coccinea) 


Lauraceae-S.E. Asia 
(Lauraceae spp.) 


Litsea (Litsea spp.) 


Rhododendron (Rhododendron . 
catawbtiense } 


Madrone (Arbutus menztesit) 


Arrowleaf balsam root 
(Balsamorhtiaa sagittata) 


Oregon grape (Berberis repans) 
Service berry (Amelanchter spp.) 
Messmate stringybark 

(Eucalyptus obliqua) 
Conifer Leaves and Needles 


Western larch (Larix vecidentalis) 
Redwood (Sequoia sempervirens) 


Western red cedar (Thuja plicata) 


Chamise (Adenostoma fasciculatun) 


Douglas-Fir (Pseudotsuga 
menziesit) 


Shortleaf pine (Pinus echinata) 


Lodgepole pine (Pinus contorta) 
Coulter pine (Pinus coultert) 
Ponderosa pine (Pinus ponderosa) 


Engelmann spruce (Picea 
engelmanntt) 


*When two values are given in this column and the one following it, the lower value of the 
absorptivity (and hence the larger value of pcL/a) are for newly fallen leaves. All 
other values are for weathered materials. Few of these absorptivities have actually been 
measured. Most are estimated from the appearance of the material. 
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CHAPTER 15 
RESIDUAL NUCLEAR RADIATION 


15-1. Gross phenomenology. 


a. Radiation emitted later than one minute after 
detonation is defined as residual. Residual radiation 
is composed of fission-product reactions and neutron 
reactions from fallout and of the neutron-induced 
gamma rays discussed in chapter 3. Fallout is the 
most significant contribution for most burst condi- 
tions of interest; it is separated into the early phase 
(which settles within 24 hr) and the delayed phase 
(which settles after 24 hr), Early fallout has military 
significance because exposure to such fallout may 
debilitate humans—they must be protected during 
those periods when the radiation intensity poses a 
threat. The design requires both dose rate and ac- 
cumulated dose (Glasstone, 1977). Neutron-induced 
soil activity (chap. 3) should be considered for burst 
elevations less than 75 ft/kt®%5. 


b. Effects of residual nuclear radiation depend on: 


- Weapon design 

- Weapon yield 

- Burst élevation 

- Range 

- Climatic conditions 
- Soil composition 


15-2. Early fallout. 


a. Air bursts. Except for the phenomenon of rain-out 
(para. 15-4), dispersion of radioactive debris from an 
air burst before it reaches the ground is usually suf- 
ficient to preclude a hazard of military significance. 

b. Surface bursts. Deposited radioactivity during 
the first day after detonation of a surface burst can 
extend for several hundred miles. Dose-rate contours 
for the environment one hour after burst over ideal 
surfaces can be prepared, as shown in figure 15-1, 
from the data presented in figures 15-2 through 
15-11. For reasonably level terrain, the values given 
in these figures should be multiplied by 0.7; rough 
and hilly terrain requires the data to be multiplied 
by a factor of 0.5 to 0.6. The effective wind is the 
average of the net wind vectors at the top and bot- 
tom of the stabilized cloud, as defined in figures 15-12 
and 15-13. The net wind vectors should be obtained 
from a plot of the wind vectors (hodograph) as a func- 
tion of altitude (fig. 15-14, for example). The 
hodograph should be plotted for equal increments of 
altitude. The net wind vectors are the vectors extend- 
ing from zero elevation to the points on the hodograph 
corresponding to the top and bottom of the stabilized 
cloud. For example, the net wind vectors and effec- 
tive wind for a typical blast.are superimposed on the 
hodograph of figure 15-14. 
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Figure 18-1. Idealized Early Fallout Dose-Rate Contour (DNA, 1972) 
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Figure 15-2. Downwind Distance and Fatlout Radiation at t = 1 Hr After Burst as a Function of Yield, 10-Knot Effective Wind (DNA, 1972) 
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Figure 15-3. Downwind Distance and Fallout Radiation at t = 1 Hr After Burst as a Function of Yield, 20-Knot Effective Wind (DNA, 1972) 
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Figure 15-4. Downwind Distance and Fallout Radiation att = 1 Hr After Burst as a Function of Yield, 40-Knot Effective Wind (DNA, 1972) 
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@ Figure 15-5. Maximum Width and Fallout Radiation att = 1 Hr After Burst as a Function of Yield, 10-Knot Effective Wind (DNA, 1972) 
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Figure 15-6. Maximum Width and Fallout Radiation at t = 1 Hr After Burst as a Function of Yield, 20-Knot Effective Wind (DNA, 1972) | 
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& Figure 15-7, Maximum Width and Fallout Radiation at t = 1 Hr After Burst as a Function of Yield, 40-Knot Effective Wind (DNA, 
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Figure 15-8. Distance to Maximum Width and Fallout Radiation at t = 1 Hr After Burst as a Function of Yield, 10-Knot Effective Wind 
(DNA, 1972) 
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Figure 15-9. Distance to Maximum Width and Fallout Radiation at t = 1 Hr After Burst as a Function of Yield, 20-Knot Effective Wind 
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Figure 15-10. Distance to Maximum Width and Fallout Radiation at t = 1 Hr After Burst as a Function of Yield, 40-Knot Effective Wind 
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Figure 15-11. Ground Zero Width and Fallout Radiation at t= 1 Hr After Burst as a Function of Yield (DNA, 1972) 
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Figure 15-14. Typical Hodograph of a Vertical Wind Structure 


c. Transition height bursts. For burst heights be- 
tween the surface and the scaled height 100 ft/kt/3 
(30.5 m/kt!’3), fallout contamination will be a 
variable problem. The elevations below which fallout 
should be considered are presented in figure 15-15. 
An estimate of the dose rate should be obtained by 
multiplying the data in figures 15~2 through 15-11 
by the adjustment factors presented in figure 15-16. 
For bursts in the upper one-quarter of the transition 
zone, the neutron-induced soil activity described in 
paragraph 3-3/represents a significant contribution 
to the residual radiation environment that must be 
considered in addition to the fallout. 


1§-3. Late fallout 


a. Decay rate. The decay rate of radioactivity nor- 
malized to the quantitative data of radiation dose 
rates presented in figures 15-2 through 15-11 is sum 
marized in figure 15-17 for any time after detonation. 
The total dosage that will be accumulated for ex- 
posure starting one hour after burst can be calculated 
from figure 15-18 and the one-hour dose rate data of 
figures 15-2 through 15-11. 

b. Mean decay rate. For practical purposes, the 
fallout radiation can be assumed to obey the t-)-? law 
in which the dosage rates at t = 1 hour are modified 


15-14 


by the factor t~!2 to obtain the dosage rates at an ar- 
bitrary time. 


15-4. Precipitation effects. 


In inclement weather, air bursts (which usually: pro- 
duce residual radiation effects of negligible military 
significance) may cause dangerous contamination on 
the ground. The effect is called precipitation scaveng- 
ing, or rain-out, and its importance depends on 
whether the burst occurs before or during precipita- 
tion; on the weapon yield, size and shape of the 
nuclear cloud, and radioactive decay rate; on the rain 
cloud size and type and duration of precipitation; and 
on the effects of precipitation on the contamination 
after it reaches the ground. Specific details regarding 
precipitation effects should be obtained from the DNA 
(1972) report. 


15-5. Uncertainty. 


a. The data presented can be considered only a guide 
for estimating residual contamination. The greatest 
uncertainty originates with the weather itself. 
Therefore, the degree of surface residual radioactive 
contamination is largely dependent on precipitation 
and wind conditions during and after detonation. 
Other effects that should be considered are weapon 
yield and HOB. 
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b, The effect of yield and HOB uncertainties on pro- 
ducing the dose rate uncertainty can be computed 
from 


1 
2. 2 
Q?, = 92 + a ; 
°F ye’ 
pn + B(log 2 - y)2/02 + p2Q2 
wr Yew vi 
(15-1) 
where 
1824, W < 9.2 R296 kt 
B = 
0.824, W > 9.2 R29 kt 
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(0.273, Z/W0-35 > 337 ft/kt0-35 


B =< 

0: Z/W°-35 < 337 ft/kt0.35 
y = 2.11 , 
n = 0.35 


Q, and Qy, are the uncertainties of HOB and yield, 
respectively; and Q- is associated with the uncer- 
tainty in the wind. It is suggested that radiation con- 
tours similar to that shown in figure 15-1 be con- 
structed for wind vectors corresponding to the mean 
value of the winds aloft data and for the means plus 
their corresponding standard deviations. From these 
data two radiation dose rates can be determined—one 
for the mean wind condition (call it B), and another 
for the mean plus the standard deviation wind condi- 
tion (say D,). Therefore, compute 


o, = Par 5 


= (15-2) 
D 
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APPENDIX A 


GLOSSARY 


we 


we 


Airblast: 


Airblast-induced Ground Shock: 


Alpha Particle: 

Azimuthal Magnetic Field: 
Beta Particle: 

Burst: 

Burst Configuration: 


Burst Elevation: 

Compton Current: 
Continuous Coverage Region: 
Crater: 


Apparent Crater: 


True Crater: 


Crater-induced Ground Shock: . 


Cratering Efficiency: 

Debris: 

Decay: 

Delivery Accuracy: 
Deposition Region: 

Depth of Burst (DOB): 
Direct-induced Ground Shock: 
Dynamic Pressure: 

Ejecta: 


Ejecta Dust Cloud: 


Electromagnetic Pulse (EMP): © 


Engulfment Radius: 


The overpressure signal following the shock wave 
driven through the air by an explosion. 

Ground response caused by passage of the airblast 
over the ground surface. : 

A helium atom stripped of its two electrons; a helium 
nucleus. 

A magnetic field whose field lines are circles in planes 
parallel to the ground surface. 


_ A free electron, 


Explosive release of energy, explosion. 

Description of the location of the explosive center of 
growth relative to the ground plane, and to the 
target or observation point. 

Height of burst; distance from the ground surface up 
to the explosive CG. 

The current resulting from the collision of gamma 
rays with air molecules and the outward displace- 
ment of electrons detached by the collisions. 

The area where sufficient ejecta is deposited to pro- 
vide complete coverage of the ground surface. 


The open hole left after a cratering explosion, smaller 
than the true crater due to partial filling by 
fallback. 

The excavation or hole defining the limit of dissocia- 
tion (separation) of the natural site material. 

Late time stress and particle motion signals at- 
tributed to crater formation. 

The ratio of crater volume to weapon yield for crater- 
ing explosions. 

Ejecta deposited on the ground surface; wreckage 
from damaged structures or vegetation. 

Decrease in amplitude; also, radioactive decay, the 
natural fission process which leads eventually to 
the formation of stable (nonradioactive) substances. 

The CEP of a weapon system or delivery system. 

The source region for EMP; the region where the 
primary gamma ray interactions with the at- 
mosphere occur. 

The distance from the ground surface down to the ex- 
plosive center of gravity; the depth of burial. 

Ground response resulting from the coupling of 
energy in the ground in the vicinity of the crater. 

One-half the product of the density and square of the 
particle velocity of air behind the shock front. 

Material thrown out of the region surrounding ground 
zero during crater formation. 

The portion of ejecta particles forming a dust cloud 
subject to prevailing wind forces. 

The “radio flash” on broad band electric and 
associated magnetic fields resulting from nuclear 
explosions. 

The maximum radius of the fireball. 
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Enhanced Weapon: 


Equivalent Ejecta Thickness: 


Equivalent Yield: 


Fallback: 

Fallout: 

Fireball: . 
Fission: 

Fluence Spectrum: 
Free Field: 


Fuel Density: 
Fusion: 


Gamma Rays: 


Ground Shock: 
Ground-Surface Dust Cloud: 


Ground Zero (GZ): 
Height of Burst (HOB): 


“High Yield” Weapon: 
Hydrodynamic: 


Impedance : 
Impulse: 


Laydown Pattern: ~ 


Load: 


A-2 


A weapon designed to enhance some aspect of bomb 
performance, e.g., the neutron bomb or an earth 
penetrating reentry vehicle/warhead, which respec- 
tively enhance neutron radiation at the expense of 
blast and ground shock and ground shock at the ex- 
pense of airblast and radiation. 

An ejecta depth based on uniform azimuthal distribu- 
tion which ignores the normally observed ray pat- 
tern and azimuthal variation. 

The yield of a weapon detonated under different con- 
ditions that would produce the same effect at the 
point of observation. 

Material lofted by the cratering process which does 
not escape but drops back into the crater. ; 
Radioactive bomb case debris and dust which settles 

out of the dust cloud. 

The high temperature, high pressure expanding 
plasma core produced by a nuclear explosion. 

The splitting of unstable heavy element atoms accom- 
panied by the release of energy. 

A spectrum relating radiation intensity (fluence) to 
energy level. 

Locations in the air or ground unaffected by the 
presence of man-made structures. 

A measure of the availability of fuel for evaluating 
firestorm potential, i.e., weight of fuel per unit 
ground surface area. 

The combination of light nuclei into heavier nuclei 
accompanied by the release of energy; the basis of 
thermonuclear (hydrogen) bombs. 

Radiation similar to visible light or x-rays but with 
very short wavelengths; hard x-rays. 

The ground stress and ground motion resulting from 
direct energy deposition and/or airblast’ loading 
from an explosion. 

The portion of the low-lying dust cloud comprising 
dust that originates from the ground surface (out- 
side the crater). 

The point on the ground surface directly above or 
below the burst point. 

The elevation of the explosion’s center of gravity 
above the ground surface. 

A general class of nuclear weapons that release a 
major portion of their energy as radiation; high 
yield to mass ratio devices. 

Fluid like; pertaining to the expansion of the hot 
plasma core during the aed phases of a nuclear 
explosion. 

The product of material density, e, and a 
characteristic wave speed, c. 

The integral of the overpressure, usually for the 
positive phase duration. 

An attack strategy utilizing the delivery of a number 
of weapons for subsequent simultaneous detona- 
tion. Delivery may be in an array to attack point, 
line, or areal targets. 

The integral of the product of the ‘ial concentration 
and the air-intake-volume rate over’a lapsed time. 


“ 
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Low-lying Dust Cloud: 


“Low Yield’ Weapons: 
Mach Reflection Region: 
Mach Stem: 


Million Electron Volts (MeV): 
Multiple Burst Attack: 
._Laydown Pattern: 


Sequential Attack: 


“Nail Driving:” 


Neutrons: 


Neutron-induced Gamma Activity: 


Nuclear Radiation: 


Nuclear Weapons: 
Outrunning: 


Overpressure: 
Pattern Bursts: 
Plastic Zone: 


Polar Electric Field: 


Positive Phase (duration): 
Precipitation Scavenging: 


Precursor (shock): 
Rad: 


Radial Electric Field: 


Radiant Exposure: 
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The dust cloud that remains in contact with the 
ground surface. 

A general class of weapons which release a major por- 
tion of their energy as kinetic energy of bomb case 
debris; low yield to mass ratio devices. 

The ground surface region below a low-altitude burst 
where reflected shock waves overtake and merge 
with the incident shock waves to form a Mach stem. 

The nearly vertical shock front formed by the merg- 
ing of reflected and incident shock waves in the 
Mach reflection region. 

Temperature or energy unit for high energy particles. 


An attack using simultaneous detonation of an array 
of previously delivered weapons. 

Repeated attack of a target, especially when crater- 
ing is a primary kill mechanism, or can be used to 
dig out a buried target; also known as Nail Driv- 
ing attack. 

A sequential multiple-burst attack designed to dig out 
and kill a buried target. 

A particle of neutral electrical charge found in most 
atomic nuclei; atomic weight ~ 1. 

Residual radioactivity induced in soil by neutron bom- 
bardment during a nuclear explosion. 

The alpha, beta, and gamma radiation resulting from 
a nuclear explosion. Commonly divided into intitial 
or prompt radiation (0 < t < 1 min) and residual 
radiation (t > 1 min). 

Weapons whose explosive energy is the result of fis- 
sion or fusion processes. Includes both atomic (fis- 
sion) and thermonuclear (fission/fusion) weapons. 

The arrival of ground signals propagated through 
underlying higher velocity layers at the ground sur- 
face before the arrival of the airblast shock front. 

The airblast pressure (in excess of the ambient 
pressure) resulting from an explosion. 

Multiple burst attacks utilizing simultaneous detona- 
tion of an array of weapons. 

The zone of soil/rock surrounding a crater in which 
significant plastic deformation has occurred. 

An electric field normal to the earth’s surface; a part 
of the EMP signal associated with an azimuthal 
magnetic field. 

The length of time the airblast signal shows a positive 
overpressure. 

Removal of radioactive dust from the dust cloud of a 
nuclear explosion by rain; rainout. 

A shock wave running ahead of the normal shock 
front, caused by radiative heating of thermally 
nonideal surfaces. 

A unit of absorbed radiation dose equivalent to 100 
ergs of radiation per gram of absorbing material or 
tissue. 

An electric field directed radially from the center (of 
a nuclear explosion). 

The amount of thermal radiation received per unit 


area (cal/em?). 
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Radiated Pulse: 


Radiation Dose: 
Rainout: 


Reflectance: 


Regular Reflection Region: 


Rupture Zone: 


Secondary Gamma Radiation: 


Sequential Detonation: 


Shake: 

Shock Wave: 

Slant Range: 
Spontaneous Ignition: 

- Superseismic: 
Temperature Inversion: 
Thermal Efficiency: 
Thermal Radiation: 


Thermally Ideal Surfaces: 


Thermally Nonideal Surfaces: 


Transition Height: 


Transmittance: 


- ‘Triple Point: 


Waveform: 
Weapon Design: 


X-ray: 


Yield: 


That portion of the EMP disturbance radiated by the 
induced electric currents. 

Total quantity of ionizing (nuclear) radiation. 

Removal of radioactive dust from a nuclear dust cloud 
by precipitation; precipitation scavenging; rain- 
induced fallout. 

Ratio of the radiant energy reaching a surface facing 
the fireball to that which would reach the same sur- 
face in the absence of clouds, the ground, and other 
reflective surfaces. 

The region on the ground surface below a low altitude 
burst where the reflected wave does not overtake 
and merge with the incident shock; the region 
within a circle of radius equal to the height of burst, 
centered at ground zero. 

A region surrounding the true crater in which crack- 
ing and crushing occur. It, in turn, is surrounded 
by the plastic zone. 

Gamma radiation produced by neutron interaction 
with materials outside the nuclear weapon. 

An attack strategy in which a target is repeatedly at- 
tacked, such as the “nail driving” attack to dig out 
and kill a. buried target. 

10-8 second 

The air-blast wave propagated by an explosion, 
characterized by discontinuous rise to peak 
pressure. 

Spherical range from the explosive center of gravity 
to the point of interest. 

Ignition as the result of radiant exposure to the ther- 
mal pulse from a nuclear explosion. 

Ground signals moving faster than the airblast shock 
front. 

Layer of air in which the wemnpereture increases with 
height. 

Ratio of the thermal energy deaieted by the fireball 
to the total energy release or yield of the weapon. 

Radiation from the fireball in the ultra-violet, visi- 
ble, infrared spectral bands. 

Perfectly reflective, nonabsorbing surface, e.g., water, 
ice. 

Absorbtive, nonreflective surface. 

The minimum burst elevation for which dust pickup 
and fallout do not occur. 

Ratio of the radiant energy reaching a surface facing 
the fireball to that which would reach the surface 
through a perfectly transparent atmosphere. 

The point at which the Mach stem and incident 
reflected shocks intersect. 

The time-dependent amplitude history of the signal. 

Weapon construction factors influencing the perfor- 
‘mance of a nuclear weapon; cf’ enhanced weapon. 

A gamma_-like radiation not associated with nuclear 
explosions. 

Explosive energy released ee a nuclear explosion, 
commonly expressed in terms of the weight of an 
energy equivalent TNT charge. 
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